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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 

Requirements  for  the  Degree  of  Doctor  of  Philosophy 

LIGHT  SCATTERING  STUDIES  OF 
CLATHRATE  HYDRATE  FORMING  SYSTEMS 

By 

Ralph  Karl  Zachary 
May,  1998 

Chairperson:  Dr.  Samuel  0.  Colgate 
Major  Department:  Chemistry 

The  current  importance  of  clathrate  hydrates  lies 
primarily  in  the  problems  they  pose  for  the  oil  and  gas 
industry;  in  the  future  they  may  represent  an  important  source 
of  energy.  These  studies  were  persued  in  an  effort  to 
determine  the  utility  of  light  scattering  for  elucidating 
changes  in  the  dynamical  structure  of  aqueous  solutions  that 
occur  prior  to  hydrate  formation.  Such  an  understanding  has 
direct  relevance  to  the  kinetics  of  clathrate  hydrate 
formation,  and  may  contribute  to  our  knowledge  of  the  many 
other  systems  in  which  water  plays  an  important  role. 

Aqueous  solutions  of  tetrahydrofuran  were  at  the  center 
of  this  research  which  included  experimental  spectroscopic 
studies  of  such  solutions  as  a function  of  temperature  and 
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system  composition.  Light  scattering  spectra  were  found  to 
contain  important  information  concerning  the  extent  and  nature 
of  hydrogen  bonding,  and  the  structural  changes  that  precede 
hydrate  nucleation.  Specifically,  the  OH  stretching  band  due 
to  water  and  the  COC  stretch  band  due  to  tetrahydrofuran  were 
seen  to  vary  in  a predictable  manner  with  changes  in  the 
composition,  temperature,  and  thermal  history  of  the  system. 

Computer  simulation  studies  were  performed  and  used  to 
aid  in  the  interpretation  of  experimental  spectra.  Low 
waveniimber  polarizability  fluctuations  were  calculated  and 
found  to  give  reorientation  times  in  good  agreement  with  those 
found  using  other  methods. 
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CHAPTER  1 

THE  PHYSICAL  CHEMISTRY  OF  CLATHRATE  HYDRATES 


Their  Occurrence  and  Significance 

A wide  variety  of  chemical  species  interact  with  water  at 
elevated  pressures  to  form  a class  of  solids  known  as 
clathrate  hydrates.  Though  similar  in  appearance  to  ice,  they 
are  readily  distinguished  by  their  chemical  and  physical 
properties,  particularly  in  their  appearance  well  above  the 
normal  freezing  point  of  water.  The  fundamental  origin  of 
these  differences  lies  in  the  fact  that  hydrates  possess  two 
distinct  kinds  of  chemical  interactions:  the  ice-like  water- 
water  interactions  of  the  lattice  framework,  dominated  by 
hydrogen  bonding,  and  water-guest  interactions  which  are  due 
almost  entirely  to  van  der  Waals  interactions^. 

Although  this  unusual  combination  of  interactions  makes 
hydrates  interesting  from  a theoretical  standpoint,  current 
interest  in  their  properties  is  primarily  the  result  of  their 
two-fold  practical  significance.  First,  several  components  of 
natural  gas  make  excellent  hydrate  guests  and  the  formation  of 
solid  hydrates  in  natural  gas  pipelines  escalates 
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transportation  costs  dramatically.  Second,  evidence  is 
rapidly  accumulating  to  suggest  that  naturally  occuring 
hydrates  represent  a potentially  enormous  resevoir  of  natural 
gas^.  In  fact,  recent  geological  surveys  suggest  that  the 
total  amount  of  oxidizable  carbon  trapped  in  naturally 
occuring  hydrates  may  surpass  the  energy  content  of  the  total 
conventional  fossil  fuel  reserves  by  as  much  as  a factor  of 
two^"'* . 

The  term  "clathrate"  was  coined  by  Powell  for  "the 
structural  combination  of  two  substances  which  remain 
associated  not  through  strong  attractive  forces  but  because 
strong  mutual  binding  of  the  molecules  of  one  sort  makes 
possible  the  firm  enclosure  of  the  other. In  clathrate 
hydrates  the  lattice  is  composed  of  water  molecules  and  they 
are  bound  to  one  another  by  hydrogen  bonding. 

Clathrate  hydrates  are  thus  distinct  from  other  types  of 
hydrates  in  that  they  are  of  the  host-guest  type.  They  may  be 
distinguished  further  by  the  weakness  of  the  interaction 
between  guest  and  host,  as  opposed  to  the  semi-clathrates 
which  contain  significant  hydrogen-bonding  interactions 
between  guest  and  host,  and  as  opposed  to  ionic  hydrates  with 
their  large  electrostatic  interactions®.  These  distinctions 
are  represented  schematically  in  Figure  1-1. 
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Figure  1-1.  Schematic  classification  of  hydrates. 
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Classification  of  Hydrates 

Clathrate  hydrates  are  classed  among  inclusion  compounds, 
along  with  lemellar  and  channel  structures.  Structurally  they 
consist  of  guest  molecules  trapped  in  cages  of  high  symmetry 
which  are  formed  by  the  host  (water)  lattice. 

Thermodynamically  they  may  be  viewed  as  solid  solutions  of 
guests  within  a host  lattice.  Although  they  are  thus 
necessarily  nonstoichiometric,  the  variation  in  guest 
stoichiometry  is  restricted  to  a range  limited  at  one  end  by 
the  minimum  cage  occupancy  required  for  thermodynamic 
stability  with  respect  to  ice  or  aqueous  solutions  and,  at  the 
other  end,  by  complete  filling  of  the  appropriately  sized 
cages.  This  range  is  significant  but  rather  narrow  and 
complicated  early  efforts  to  characterize  hydrates^. 

Clathrate  hydrates  have  been  classified  according  to  the 
chemical  nature  of  the  guest  species,  by  their  thermodynamic 
properties,  and  by  their  crystalline  structure.  The  weakness 
of  the  guest-host  interaction  suggests  that  the  fundamental 
classification  is  according  to  the  structure  of  the  water 
lattice  and  so  this  approach  will  be  described  first. 
Classification  bv  Clathrate  Structure 

With  only  two  exceptions,  all  known  clathrate  hydrates 
have  water-lattices  corresponding  to  one  of  three  types 
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denoted  Structure  I (SI),  Structure  II  (SII),  and  the  recently 
discovered  Structure  H (SH)®'^.  Their  geometry  is 

characterized  by  specific  combinations  of  nearly  spherical 
voids.  There  are  five  types  of  these  voids  and  all  three 
structures  contain  pentagonal  dodecahedra  which,  because  they 
consist  of  12  pentagonal  faces,  are  frequently  denoted  5^^. 
Structure  I hydrates  also  contain  tetrakaidecahedra,  denoted 
6^5’'-^  to  indicate  their  construction  from  two  hexagonal  and 
twelve  pentagonal  faces.  The  second  cage  type  in  Structure  II 
is  the  hexakaidecahedron  denoted  6^5^^.  Structure  H consists 
of  three  cage  types:  pentagonal  dodecahedra,  the  20-hedron 
denoted  6®5^^,  and  the  12-hedron  4®5®6®.  The  five  cage  types 
are  illustrated  in  Figure  1-2  and  their  geometric  properties 
summarized  in  Table  1-1. 

The  water  lattice  of  Structure  I hydrates  is  body- 
centered  cubic  and  belongs  to  the  space  group  Pm3n.  The  unit 
cell  dimension  is  typically  12  A and  does  not  vary 
significantly  with  variations  in  the  guest  molecule.  A single 
unit  cell  has  the  formula  2Mi»6M2»  46  H2O,  where  Mi  and  M2  refer 
to  the  5^^  and  cavities,  respectively. 

Structure  II  hydrates  posses  a diamond  cubic  water 
lattice  that  belongs  to  the  space  group  Fd3m.  Typical  cell 

o 

dimensions  are  17  A on  a side  and,  again,  do  not  depend 
strongly  on  the  nature  of  the  guest.  The  unit  cell  has  the 
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formula  16Mi»8M3*136  H2O/  with  and  M3  the  5^^  and  cageS/ 

respectively. 


Table  1-1.  Cavities  occuring  in  clathrate  hydrates 


Cavity  Type 

5^2 

5^262 

5^254 

435663 

5^268 

Average  Cavity 
Radius,  A 

3.90 

4.33 

4.68 

4.06 

5.71 

Coordination  Number 

20 

24 

28 

20 

36 

Cavities/  Unit  Cell 
(SI) 

2 

6 

— — — 

— 

— 

Cavities/  Unit  Cell 
(SII) 

16 

— 

8 

— 

— 

Cavities/  Unit  Cell 
(SH) 

3 

— 

— 

2 

1 

Table  1-2.  Properties  of  the  Water  Lattices  in  Clathrates 


Structure 

I 

II 

Crystal  System 

Cubic 

Cubic 

Space  Group 

Pm3n  (#223) 

Fd3m  (#227) 

Lattice  Description 

Body  Centered 

Diamond 

Lattice  Parameter, 

12 

17.3 

Ideal  Unit  Cell 
Formula 

2Mi»6M2«  4 6 H2O 

16Mi»8M3»13  6 H2O 
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Figure  1-3 (a) . The  SI  unit  cell. 


Figure  1-3 (b) . Projection  of  the 
SI  lattice. 
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Figure  1-3 (a) . 


The  SII  unit  cell. 


Figure  1-3 (b)  Alternative  view  of 
the  SII  unit  cell. 
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Figure  1-5.  The  Structure  H unit  cell. 
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structure  H hydrates  have  only  recently  been  discovered 
and  contain  a water  lattice  which  is  hexagonal  in  symmetry. 
The  unit  cell  has  the  formula  3Mi*2M4»M5»34H20  where  M^,  M4,  and 
M5  are  the  5^^,  and  cavities  respectively. 
Structure  H hydrates  are  of  only  incidental  interest  to  this 
study. 

The  unit  cells  of  the  three  hydrate  types  are  illustrated 
in  figures  1-3,  1-4,  and  1-5.  Several  of  their  more  important 
features  are  summarized  in  Table  1-2. 

Clathrate  hydrates  resemble  ice  to  the  naked  eye  and  do 
so  at  the  molecular  level  in  the  sense  that  each  water 
molecule  is  hydrogen-bonded  to  four  nearest  neighbors,  as  in 
all  the  known  forms  of  ice.  The  guest  molecules  play  a 
stabilizing  and  space-filling  role  which  does  not  depend 
strongly  on  their  identity.  The  average  departure  of  0-0-0 
angles  from  the  tetrahedral  values  found  in  ices  and  is 
only  3.7°  and  3.0°  in  structures  I and  II  respectively.  This 
is  in  contrast  to  the  rather  large  deviations  found  in  the 
high-pressure  ices  II,  III,  V and  VI.  The  0-0  bond  lengths  on 
average  exceed  those  in  hexagonal  ice,  where  they  are  2.76  A, 
by  only  1%  and  are  comparable  to  those  in  most  of  the  other 
forms  of  ice^°"^^. 

Aside  from  the  cages,  the  most  obvious  structural  feature 
distinguishing  clathrates  from  the  various  forms  of  ice  is  the 
predominance  of  planar  (structure  II)  or  nearly  planar 
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(structure  I)  five-membered  rings.  The  less  numerous 
pentagons  which  occur  in  ice  IX  and  III  are  far  from  planar. 
The  hexagonal  rings  are  also  planar  in  structure  I or  very 
nearly  so  in  the  case  of  structure  II.  This  is  in  marked 
contrast  to  the  puckered  hexagons  found  in  other  forms  of  ice. 
This  planarity  means  that  the  three  water  molecules  hydrogen- 
bonded  to  the  molecule  at  each  end  of  a hydrogen  bond  are 
eclipsed  when  viewed  along  that  bond,  or,  to  use  Bjerrum's 
terminology,  that  the  arrangement  about  the  bond  is  mirror 
symmetric  (ms)  . In  cubic  ice  Ic,  and  probably  in  ices  VII  and 
VIII,  the  orientations  about  all  bonds  are  staggered  or  trans 
in  a center  symmetric  (cs)  arrangement.  This  is  the  case  in 
ice  Ih,  except  for  one-fourth  of  the  bonds,  those  which  lie 
parallel  to  the  c axis,  and  about  which  the  arrangement  is 
also  . 

Bjerrum's  analysis  hypothesized  that  the  ms  arrangement 
was  more  stable  than  the  cs  but  that  it  was  probably  not 
possible  to  construct  a crystal  in  which  more  than  one-fourth 
of  the  bonds  were  of  the  ms  type.  The  requirements  of  ring 
and  cage  closure  produce  some  distortions  of  most  bonds  in 
clathrate  hydrates  from  the  ideal  tetrahedral  ms  arrangements 
found  around  c bonds  in  Ih.  This  distortion  of  angles  and 
bond  lengths  and  the  weaker  attractive  forces  between 
nonbonded  neighbors,  rather  than  any  fundamental  instability 
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of  ms  with  respect  to  cs  bonding,  is  what  accounts  for  the  low 
stability  of  the  empty  clathrate  hydrate  lattices  with  respect 
to 

Pentagonal  rings  outniamber  hexagonal  rings  in  a ratio  of 
8:1  in  structure  I and  9:1  in  structure  II.  Both  are  close  to 
the  ratio  of  8.62:1  which  would  make  the  average  angle 
tetrahedral  if  all  faces  were  exactly  planar.  The  departures 
from  planarity  which  do  occur  are  in  the  direction  of 
equalization  of  hydrogen  bond  lengths  throughout  the 
structures . 

Thermodynamic  Classification 

A second  scheme  for  the  classification  of  clathrate 
hydrates  is  based  on  the  phase  behavior  of  hydrate-forming 
systems.  It  is  a convenient  way  of  describing  the  conditions 
under  which  a hydrate  is  stable.  Most  hydrate  forming  systems 
have  a phase  diagram  similar  to  that  of  SO2/H2O  illustrated  in 
figure  1-6  in  which  h stands  for  hydrate,  g for  gas,  I for 
ice,  li  for  the  liquid  phase  rich  in  water  and  I2  the  liquid 
phase  rich  in  hydrate  former.  Systems  with  such  a diagram  are 
denoted  i2.  The  region  of  hydrate  stability  is  bounded  by  the 
h-I-g,  h-l^-g  and  h-li-l2  equilibrium  lines^'*. 

Two  invariant  points  Qi  and  Q2  mark  the  coexistence  of  h, 
I,  li,  and  g on  the  one  hand  and  of  h,  1^,  I2  and  g on  the 
other.  This  phase  behavior  reflects  the  sparing  solubility  of 
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most  hydrate-formers  in  water  and  such  systems  are 
characterized  by  the  temperatures  and  total  pressures  at  the 
quadruple  points  Qi  and  Q2. 

Among  the  guests  which  form  i2  systems  are  the  SI  formers 
CI2,  CO2,  N2O,  H2S,  SO2,  PH3,  ASH3,  C2H2,  C2H6,  and  CH3CI.  SII 

formers  include  SFg,  C3Hg,  iso-butane,  cyclopentene, 
cyclopentane,  and  CH2CI2. 

The  Q2  point  is  absent  for  light  gases  whose  h-l2-g 
equilibrium  line  ends  at  the  l2~q  critical  point  in  the 
presence  of  hydrate  before  it  intersects  the  h-l^-g.  Such 
systems  are  denoted  ii,  and  in  them  hydrate  may  exist  to 
indefinitely  high  temperatures  so  long  as  a pressure  is 
maintained  that  is  at  least  equal  to  the  steeply  rising 
decomposition  pressure. 

The  most  important  species  forming  i^  systems  with  water 
all  form  SI  hydrates  and  include  the  noble  gases  Ar,  Kr,  and 
Xe  in  addition  to  N2,  O2,  CH3,  and  C2H4. 

If  the  guest  species  is  sufficiently  soluble  in  water 
there  is  only  one  liquid  phase  and  the  quadruple  point  Q2  is 
absent.  For  some  hydrates  it  is  replaced  (as  the  maximum 
temperature  of  coexistence  of  hydrate,  liquid,  and  gas)  by  the 
congruent  melting  point.  The  phase  relationships  in  such  a 
case  are  more  conveniently  shown  in  the  temperature- 
composition  projection  of  the  phase  diagram.  In  figure  1-7 
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such  a project  is  shown  for  the  SI  forming  ethylene-oxide 
water  system^^. 

Guest  Species  of  SI  and  SII  Clathrate  Hydrates 

Some  100  chemical  species  are  known  to  form  clathrate 
hydrates.  Suitability  of  a guest  is  determined  primarily  by 
its  size:  its  molecules  must  be  small  enough  to  be  contained 
by  the  cages  available  in  a particular  hydrate  structure.  The 
chemical  identity  of  the  species  is  relevant  only  in  so  far  as 
the  molecules  must  not  contain  moieties  leading  to  covalent  or 
ionic  bonding  with  water. 

While  each  cage  of  a hydrate  may  contain  no  more  than  one 
guest,  they  need  not  necessarily  all  contain  the  same  kind  of 
guest.  This  leads  to  the  distinction  between  simple 
hydrates,  those  containing  only  one  type  of  guest,  and  the  more 
complex  double-  and  mixed-hydrates^®.  In  the  latter,  the 
various  guest  species  occupy  cages  of  all  available  types 
while  in  the  former  guest  species  of  a particular  type  are 
found  only  in  one  type  of  cage. 

The  Dynamics  of  Hydrates 

In  the  years  which  immediately  followed  the  1811 
discovery  of  hydrates  by  Sir  Humphrey  Davy^^,  most  research  in 
the  field  was  devoted  to  the  determination  of  the  composition 


15 


of  various  hydrates.  This  was  followed  by  a period  of  intense 
interest  in  the  structure  characteristics  of  clathrates  and 
the  phase  behavior  of  clathrate-forming  systems^®.  Research 
along  these  lines  lead  to  the  formulation  of  an  elegant 
statistical  mechanical  model  of  hydrates  by  van  der  Waals  and 
Platteeuw^®'^°;  a model  that  has  been  particularly  useful  for 
the  oil  and  gas  industry  in  their  efforts  to  predict  the 
thermodynamic  conditions  conducive  to  hydrate  formation. 
Recently  interest  has  shifted  toward  a more  detailed  study  of 
the  dynamics  of  hydrates,  and  into  the  kinetics  of  their 
formation. 

Lattice  Dynamics 

In  common  with  hexagonal  ice,  and  the  other  disordered 
forms  of  ice,  the  orientational  disorder  of  the  water 
molecules  in  clathrate  hydrates  manifests  itself  in  a 
rotational  mobility  affecting  the  dielectric  and  nuclear 
magnetic  resonance  properties  at  relatively  high 
temperatures^^.  Reorientations  of  the  water  molecules  result 
in  interchanges  among  the  ~(3/2)N  lattice  configurations  until 
the  temperature  becomes  so  low  that  the  rate  of  interchange  is 
negligible  and  the  configurational  disorder, 
corresponding  residual  entropy,  is  "frozen-in." 


with  its 
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Figure  1-7.  Typical  phase  diagram  for  a system  containing  a 
water  soluble  hydrate  former. 
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Dielectric  measurements,  in  particular,  suggest 
significant  dynamical  differences  between  the  ices  and 
clathrate  hydrates.  The  Cole-Cole  complex  permittivity  locus 
of  ice  Ih is  accurately  a semi-circle;  i.e.,  reorientation  of 
the  water  molecules  occurs  with  a single  relaxation  time^^. 
The  temperature  dependence  of  this  reorientation  time  follows 
Arhenius-like  behavior  with  an  activation  energy  of 
approximately  13.2  kcal  mol'^.  The  classic  analysis  of  the 
dielectric  properties  of  ice  by  Bjerrum  models  the 
reorientation  of  water  molecules  in  the  ice  lattice  as  the 
diffusion  of  defects  originating  instrisically  from  thermal 
excitation^^. 

In  contrast,  the  structure  II  hydrates  have  a complex 
permitivity  locus  described  in  terms  of  a spread  of  relaxation 

times  about  Tq,  the  most  probable  relaxation  time,  described 
quantitatively  by  the  parameter  a in 

e*  = g'  - iq"  = G„^  + (Gq  - G^J  / [1  + (icoTp)  [1-1] 

In  the  case  of  ice  Ih,  a = 0. 

Furthermore,  denoting  the  thermodynamic  temperature  by  T, 

the  1/T  dependence  of  log  Xq  deviates  significantly  from 

linearity  which  is  probably  due  to  the  increasing  predominance 
at  low  temperatures  of  relaxation  controlled  by  Bjerrum 
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defects  of  extrinsic  origin.  Although  there  appears  to  be  a 

loose  correlation  of  Tq  with  guest  polarity,  the  correlation 

is  not  close,  suggesting  that  these  differences  must  be,  at 
least  in  part,  attributed  to  differences  between  the  water 
lattices  in  ice  and  in  clathrates^'’ . 

Results  from  both  proton  nuclear  magnetic  resonance  (NMR) 
spectra  and  dielectric  measurements  indicate  a second 
characteristic  distinguishing  the  dynamics  of  water  molecules 
in  ice  from  their  motion  in  hydrates^^.  In  addition  to  the 
multiplicity  of  relaxation  processes  present  in  hydrate 
lattices  and  discussed  above,  the  rates  of  water  reorientation 
differ  greatly  in  the  two  environments.  Though  hydrate  water 
molecules  reorient  rapidly  relative  to  molecules  in  other 
solids,  those  reorientation  rates  are  only  one-half  those  in 
ice . 

This  difference  is  particularly  significant  in  that  it 
changes  the  dynamic  relationship  between  rotation  and 
translation  in  hydrates.  Water  molecules  in  ice  diffuse 
almost  an  order  of  magnitude  faster  than  they  reorient  about 
a fixed  position  in  the  crystal  structure.  Water  molecules  in 
the  hydrate  lattice  reorient  20  times  faster  than  they 
diffuse^®. 

What  remains  in  dispute  is  the  magnitude  and  kind  of 
influence  exerted  by  the  guest  molecules  over  the  rotational 
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and  translational  motion  of  water  in  hydrate  lattices . Proton 
NMR  studies  aimed  at  elucidating  the  motion  of  the  guests  do 
not  yield  evidence  of  hydrogen  bonding,  while  dielectric 
studies  seem  to  imply  the  occasional  formation  of  hydrogen 
bonds  between  guest  and  host  molecules  so  as  to  inject  Bjerrum 
defects  of  a specific  type  into  the  lattice. 

These  two  methods  seem  incapable  of  resolving  the  dispute 
since  the  extent  of  hydrogen  bonding  would  presumably  have  to 
exceed  several  percent  for  it  to  be  readily  detected  by  proton 
NMR,  while  the  dielectric  relaxation  behavior  can  be  accounted 
for  by  supposing  only  a small  fraction  of  guests  to  be 
hydrogen  bonded  (~2  X 10"®)^^. 

Guest  Dynamics 

What  is  not  in  dispute,  and  which  is  in  distinct  contrast 
to  other  solid  systems,  is  the  remarkable  degree  of  rotational 
mobility  of  the  guest  species.  The  permittivity  measured  on 
the  high-frequency  side  of  the  dispersion  region  associated 
with  reorientation  of  the  water  molecules  reaches  a limiting 
value  e„i  which  depends  on  the  nature  of  the  guest  molecule  and 
on  the  temperature.  This  limiting  value,  as  measured  for 
structure  II  hydrates  at  a fixed,  relatively  high  temperature, 
is  very  nearly  proportional  to  the  square  of  the  guest 
molecule's  dipole  moment.  Above  about  150  K the  value  varies 
linearly  with  1/T  for  most  hydrates.  Taken  together  these 
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observations  indicate  that  the  guest  molecules  are  undergoing 
reorientation  at  rates  considerably  in  excess  of  the 
frequencies  of  the  measurements,  ~1  MHz,  either  isotropically 
or  between  two  or  more  preferred  orientations  of  similar 
energy. 

A study  of  the  dielectric  absorption  associated  with 
reorientation  of  the  guest  molecules  in  the  hydrates  of 
ethylene  oxide,  tetrahydrofuran,  and  acetone  gave  estimated 
relaxation  times  on  the  order  of  about  10'^^  sec  at  90  K.  The 
Eyring  activation  energies  (E^  - RT)  were  estimated  at  0.5 
kcal  mol“^  for  ethylene-  oxide  hydrate  and  0.3  kcal  mol'^  for 
the  structure  II  hydrates. 

Also  relevant  is  the  observation  that,  within  each 
structure  type,  the  relaxation  rate  of  the  guest  decreases  and 
the  activation  energy  increases  with  increase  in  size  of  the 
guest  molecule^®. 

Hysteresis  and  the  Kinetics  of  Hydrate  Formation 

An  important  feature  of  hydrate- forming  systems,  and  one 
that  has  plagued  efforts  to  study  their  phase  behavior,  is 
the  extreme  hysteresis  and  metastability  found  in  such 
systems.  This  feature  complicates  the  precise  determination 
of  phase  boundaries  and  renders  ambiguous  the  results  from 


different  laboratories. 


H^^drate  Formation 
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Figure  1-8.  Plot  of  pressure  versus  temperature  for  a typical  hydrate  formation 
experiment . 
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Figure  1-9.  Overview  of  a proposed  hydrate  formation  process. 
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In  figure  1-8  the  hysteresis  effect  is  apparent  in  the 
isochoric  pressure  trace  of  a hydrate  formation  and  hydrate 
decomposition  cycles^®.  The  cycle  begins  with  cooling  of  a 
hydrate  cell  containing  methane  gas  and  water.  The  letter  0 
on  the  plot  indicates  the  start  of  hydrate  formation  which  was 
accompanied  by  a sharp  drop  in  the  total  pressure.  After 
hydrate  ceased  to  form  the  temperature  of  the  sample  was 
increased  until  hydrate  decomposition  is  complete,  marked  by 
an  X on  the  plot.  The  system  was  then  cooled  again  and  it  was 
found  that  with  each  successive  cooling  run  of  the  sample,  the 
temperature  at  which  hydrate  formation  begins  was  higher  than 
in  the  previous  trial.  Results  such  as  these  have  been 
offered  as  evidence  for  the  existence  of  metastable  hydrate 
precursors  at  temperatures  well  above  the  ostensible  hydrate 
decomposition  temperature. 

In  common  with  crystallization,  the  hydrate  formation 
process  may  be  divided  into  nucleation  and  growth  phases.  In 
the  nucleation  or  induction  phase,  hydrate  precursors  form  and 
disperse  within  a supersaturated  aqueous  solution.  If  these 
precursors  manage  to  grow  to  a critical  size,  the  second  phase 
begins  and  the  hydrate  nuclei  quickly  grow  into  solid 
crystals.  The  system's  history  exerts  a powerful  influence 
over  the  length  of  the  induction  time.  This  complication  is 
largely  responsible  for  the  shortage  of  reliable  quantitative 
kinetic  data  available  for  the  induction  period^®. 
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Figure  1-9  outlines  a hydrate  formation  hypothesis 
offered  by  Sloan^^  in  which  species  are  seen  to  progress  from 
water  and  small  water  clusters  [A]  through  labile  clusters  [B] 
to  metastable  agglomerates  [C]  to  stable  nuclei  [D]  which 
marks  the  end  of  the  nucleation  period  and  beginning  of  the 
rapid  growth  process. 

Heating  of  such  a system  will  drive  the  process  in 
reverse  and  the  stable  hydrate  crystals  will  dissociate.  Yet 
the  evidence  suggests  that  microscopic  species,  ranging  in 
size  from  a single  labile  cluster  [B]  to  agglomerates [C] , 
persist  in  the  solution  to  temperatures  up  to  about  20  K above 
the  temperature  at  which  bulk  hydrate  disappears. 

What  is  missing  from  this  hypothesis  is  a more  detailed 
description  of  the  clusters  responsible  for  the  metastability 
of  hydrate  forming  systems.  Research  directed  along  these 
lines  is  in  its  infancy  and  consists  almost  exclusively  of 
preliminary  NMR  and  dielectric  measurements.  The  discussion 
of  lattice  and  guest  dynamics  given  earlier  suggests  that 
these  techniques  will  necessarily  give  an  incomplete  picture 
of  the  relevant  dynamic  processes  and  need  to  be  supplemented 
by  other  methods. 

One  class  of  methods  sensitive  to  the  dynamics  of  liquids 
concerns  the  measurement  of  scattered  light.  Although  the 
related  technique  of  infrared  spectroscopy  has  been  used  in 
the  analysis  of  hydrates,  these  studies  were  restricted  to 
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solid  hydrates^^,  solid  ice^^,  pure  water^'’,  or  to  the 
equilibrium  properties  of  aqueous  solutions  not  necessarily 
relevant  to  hydrates^^. 

There  are  two  published  reports  in  the  literature  of 
Raman  scattering  from  solid  hydrates^®  and  their  results  are 
briefly  discussed  in  Chapter  5,  but  the  focus  in  those  studies 
was  on  the  pressure  dependence  of  certain  structural  features 
of  hydrates.  More  relevant  are  the  numerous  studies  of  light 
scattering  from  pure  water,  and  from  aqueous  solutions  not 
involved  in  hydrate  formation,  those  results  will  be  discussed 
in  more  detail  in  Chapter  5.  These  pages  describe  the  first 
known  attempt  to  apply  the  insights  obtained  from  such  studies 
to  hydrate- forming  aqueous  solutions. 

Specific  Guests  of  Interest 

Three  guest  molecules  figure  prominently  in  the  studies 
reported  here  and  their  chemical  structures  are  given  in 
Figure  1-10.  The  first,  methane,  is  an  important  component  of 
natural  gas  and  its  simple,  tetrahedral  geometry  makes  it 
convenient  for  modeling  hydrate  properties.  The  methane 
hydrate  has  a lattice  structure  corresponding  to  Structure  I 
and  forms  at  a temperature  of  -78.7  °C,  at  1 atmosphere 
pressure^^ . 
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Ethylene  oxide/  or  oxirane/  is  a gas  at  ordinary 
temperatures  and  pressures  and  it  also  forms  a Structure  I 
hydrate^'*.  Its  phase  diagram  is  characterized  by  a lower 
quadruple  point  at  -2.1  °C  and  1 atmosphere  pressure. 

A larger  cyclic  ether,  tetrahydrofuran,  was  chosen  for 
contrast  with  ethylene  oxide.  Tetrahydrofuran  combines 
withwater  to  form  a Structure  II  hydrate  at  4.3  °C  and  1 
atmosphere  pressure.  It  is  a liquid  at  room  temperature  and 
totally  miscible  in  water. 
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Figure  1-10.  Hydrate  guests  of  special  interest. 


CHAPTER  2 

LIGHT  SCATTERING  FROM  MOLECULAR  LIQUIDS 
Introduction:  Classical  Theory  of  Light  Scattering 

Although  a complete  description  of  the  interaction  of  light 
with  matter  requires  a quantum  mechanical  treatment,  the  modern 
understanding  of  light  scattering  can  be  said  to  date  from  Lord 
Rayleigh's  description  of  the  central  role  played  by  polarizability 
fluctuations.  With  the  development  of  kinetic-molecular  theory  a 
classical,  yet  molecular,  theory  of  light  scattering  emerged.  This 
led  to  Lommel's  prediction  of  inelastic  light  scattering,  and  the 
first  quantum  mechanical  treatment  by  Smekal.  The  effect  was 
observed  experimentally  by  C.  V.  Raman  and  K.  S.  Krishnan  in 
liquids  and,  simultaneously,  by  G.  S.  Landsberg  and  L.  I. 
Mandelstam  in  crystals  though  neither  experimental  team  was  aware 
of  the  theoretical  predictions^^. 

The  insight  provided  by  classical  theory  is  sufficiently 
valuable  to  warrant  a brief  review  of  its  essentials.  The 
discussion  will  be  limited  to  what  is  sometimes  described  as 
"spontaneous"  light  scattering  to  distinguish  it  from  resonant  and 
nonlinear  phenomena^®”^® . 
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Consider  an  individual  molecule  subject  to  an  incident  light 
wave  described  by,  E = Eq  cos  27iVot.  The  dipole  moment,  \i,  induced 
by  this  external  electromagnetic  field  is  given  by. 


Where  a is  the  molecular  polarizability  and,  in  the  general  case, 

is  a tensorial  quantity.  In  this  section  it  will  be  assumed  that 
the  electric  field  and  induced  dipole  moment  are  parallel  to  one 
another  such  that  these  three  quantities  may  be  treated  as  scalars. 

Continuing  with  the  classical  approach,  suppose  the 
polarizability  is  a phenomenological  quantity  which  depends  on  the 
distances  between  the  atomic  nuclei  at  any  given  moment.  Denoting 
by  q^  the  coordinate  describing  a particular  vibrational  mode  of  the 

molecule,  we  may  write  a = a{q^) . If  the  displacement  is  small, 

the  polarizability  can  be  expanded  in  a power  series  in  q^  around 
the  equilibrium  value  of  this  coordinate,  q^  = 0.  That  is,  we  may 
write 


= a E 


[2-1] 


[2-2] 


After  substituting  q^  = q^o  sin  (27cvt ) we  have 


a(t)  = a + a cos  2nv  t 

O V V 


[2-3] 
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The  expression  for  the  dipole  moment  becomes, 

p(t)  = Wq  Eg  cos2nVg  t + cos2nVg  t cos2nv^  t [2-4] 

or  equivalently, 

oc  E 

]i{t)  = «g  Eg  cos2nVg  t + — ^ [cos2n(Vg  - v^)  t 

^ [2-5] 

+ cos2n  (v.  + V ) t] 

The  significance  of  this  time  varying  dipole  can  also  be 
understood  within  the  classical  framework  since  classical 
electromagnetic  theory  holds  that  such  a dipole  will  emit  light  of 
the  same  frequency. 

Several  important  features  of  light  scattering  are  already 
apparent  even  in  this  minimal  treatment.  First,  the  scattered 
light  is  a result  of  forced  oscillations  of  the  molecular  dipole 
moment  induced  by  the  electromagnetic  field  of  the  incident  light 
wave.  Rayleigh  scattering,  represented  by  the  first  term  of 
equation  2-5,  which  consists  of  light  of  unaltered  frequency,  is 
due  to  the  unperturbed  polarizability  of  the  molecule. 

Furthermore,  we  may  anticipate  the  results  of  the  quantum 
mechanical  treatment,  by  inferring  the  general  features  of  Raman 
scattering  from  the  remaining  two  terms  in  equation  2-5.  The  light 
is  primarily  scattered  by  the  electronic  properties  of  the  molecule 
and  the  atomic  nuclei  are  displaced  only  slightly.  The 
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deformability  of  the  electron  cloud  in  the  imposed  electric  field 
depends  on  the  nuclear  configuration  at  any  given  moment.  As  the 
nuclei  oscillate  about  their  respective  equilibrium  positions,  the 
deformability  of  the  electron  cloud  varies  with  the  oscillation 
frequency  of  the  nuclei.  On  the  other  hand,  deformation  of  the 
electron  cloud  may  induce  vibration  of  the  nuclei. 

The  classical  treatment  correctly  predicts  that  the  intensity 
of  the  Raman  scattering  lines  is  proportional  to  the  square  of  the 
derivative  of  polarizability  with  respect  to  the  vibrational 
coordinate^®.  If  higher  terms  in  the  expansion  2-2  are  taken  into 
consideration,  the  intensities  of  overtones  and  combination 
frequencies  can  be  similarly  related  to  the  squares  of  the  higher 
derivatives . 

The  Raman  shift  is  Av  = Vq  - Vg^at  = and,  because  it  is  a 

property  of  the  molecule  independent  of  the  incident  light's 

frequency,  Vq,  it  is  treated  as  the  independent  variable.  At  this 

point  a potentially  confusing  feature  of  the  literature  should  be 
mentioned:  the  Raman  shift  is  almost  universally  reported  as  the 

wavenumber  shift,  Av  = Av/c,  but  it  is  often  denoted  simply  v or  o 
and,  even  more  unfortunately,  is  often  referred  to  as  "frequency." 
These  pages  will  reserve  <o  for  angular  frequency,  v for  frequency, 
and  V for  wavenumber. 

The  origin  of  light  scattering  may  derive  from  the 
translational  or  rotational  motion  of  molecules  to  the  much  faster 
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(and  higher  energy)  vibrational  motion  of  nuclei  within  a molecule 
and  these  modes  of  motion  are  traditionally  treated  using  different 
techniques.  What  these  techniques  share  is  a common  focus  on  the 
fluctuation  of  the  system's  polarizability.  Detailed  treatment  of 
these  different  techniques  will  be  postponed  until  after  a brief 
description  of  the  polarizability  of  a molecule. 

Quantum  Mechanical  Treatment:  the  Placzek  Approximations^Q 

Prior  to  its  being  perturbed  by  incident  light,  an  independent 
molecule  exists  in  one  of  the  stationary  states  | 1>  specified  by. 


The  introduction  of  a time  varying  electromagnetic  field,  E(t) 
= e^“°*^  + E"  induces  the  system  to  undergo  a transition  to  a 

time  dependent  state  M/i(t)  which  satisfies. 


Here,  D is  the  dipole  moment  operator.  The  time  dependent  states 
may  be  written. 


\1>  = \1> 


1 = 1,2,3... 


[2-6] 


[2-7] 


|il;^(t)>  = ( |J>  + I4);>)  + ( |J>  + !(!);>) 


[2-8] 
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where  the  states  I4)*>  are  time  independent.  A relatively  simple 

result  is  obtained  if  one  neglects  terms  higher  than  first  order  in 
E,uses  the  completeness  of  the  set  of  stationary  states,  and 

assumes  that  the  laser  frequency,  ©0/  is  not  too  close  to  the 

frequency  of  transition  ©„i: 


|m> 


<m\  [D  'E(t)  \1> 


[2-9] 


The  shift  in  frequency  characterizing  Raman  scattering  is  due  to 
transitions  between  two  quantum  stationary  states.  The  induced- 
transition  dipole  involved  in  Raman  scattering  associated  with  a 
transition  from  state  | 1>  to  state  |m>  is  obtained  from. 


[2-10] 


After  neglecting  the  part  of  the  dipole  moment  responsible  for  the 
spontaneous  emission,  we  can  see  that  the  Cartesian  component  of 
the  induced  dipole  associated  with  this  transition  is  given  by,  to 
first  order  in  the  electric  field. 


Dfu)  - E*  [<(“„)  K 


+ [a“^(©g)  ]*  £:;  ] 


[2-11] 


34 


In  equation  2-11,  a'^jk(Wo)  is  the  jk  component  of  the  polarizability 
tensor  found  from: 


= E.  ( 


<m  £1.  nXn  £)^|  J> 

' j ' ' Jcl 

+ %) 


<m  D.  nXn  £>,  1> 

■#5  (q  + CO.) 
nm  0 ' 


[2-12] 


In  equation  2-12,  the  case  where  co„,i  < 0 corresponds  to  Stokes 

scattering,  whereas  o)„i  > 0 corresponds  to  anti-Stokes  scattering. 

The  intensity  of  the  scattered  light  is  obtained  by 
substituting  the  induced  transition  dipole  moment  into  Maxwell's 
ecjuations  and  then  solving  the  resulting  inhmogenous  wave  ecjuation. 
For  incident  laser  light  polarized  along  the  k direction,  and 
scattered  light  observed  in  the  j direction  of  polarization,  the 
result  is: 


[2-13] 


in  which  Iq  is  the  incident  intensity  given  by. 


) e:e: 

J J 


[2-14] 


Although  this  result  is  general,  it  is  not  convenient  for  the 
description  of  scattering  from  molecules. 

When  the  transition  frequencies  for  electronic  motion  are 
large  compared  to  the  separation  in  vibrational  or  rotational 
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energies,  the  approximations  introduced  by  Placzek^^  are  appropriate 
and  transform  this  into  a form  suitable  for  the  study  of 
vibrational,  rotational,  and  translational  motions  for  many 
molecular  system. 

Placzek  theory  assumes  that  the  adiabatic  approximation  is 
valid,  and  that  the  molecule  is  in  its  (nondegenerate)  ground 
electronic  state  before  and  after  interaction  with  incident 
radiation. 

The  adiabatic  approximation  permits'*^  us  to  write  the  wave 
function  of  stationary  state  | 1>  in  a form  that  separates  nuclear 
from  electronic  motion: 


il/j  {x,X^r  = d)j  ix;X) 


[2-15] 


Here  x and  X refer  to  the  electronic  and  nuclear  coordinates, 
respectively;  (|)i(x;X)  is  the  electronic  wave  function  with  nuclei 

fixed  in  a configuration  specified  by  X;  and  ^ la  (X)  is  the  nuclear 

wave  function  for  nuclei  in  the  state  specified  by  the  vibrational 
and  rotational  quantiom  numbers  collectively  represented  by  a,  and 
in  which  the  electronic  state  of  the  molecule  is  1. 

The  energy  associated  with  state  | 1>  is  found  to  be 


la 


£°  + W,  = e,  (X) 

1 la  Ja  ' ' 


[2-16] 
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in  which  Si°  is  the  eigenenergy  of  the  electrons  with  nuclei  fixed 
in  the  configuration  X and  is  the  eigenenergy  associated  with 
the  wave  function  ^ (X)  . 

We  may  now  take  advantage  of  the  second  approximation,  and  set 
the  initial  state  1 1>  to  be  I ga>  and  the  final  state  |m>  to  be 
|gb>.  That  is,  we  have  explicitely  labeled  the  ground  electronic 
state  by  g and  have  chosen  a and  b to  label  sets  of  quantum  numbers 
describing  the  state  of  nuclear  motion.  Then  the  expression  for 
the  relevant  polarizability  is  given  by. 


<gh\D  \fcXfc\D^\ga>  <gb\D^\fcXfc\D  \ga> 

If  the  incident  frequency  Qq  is  far  from  any  electronic  absorption 

band  such  that  - o is  large  compared  to  the  shift  of  the 

electronic  levels  resulting  from  nuclear  motion  then  we  may  neglect 
©ca  ®cb  iri  the  denominator  of  equation  2-17. 

The  result  of  the  Placzek  approximations  may  be  written, 

= <b|o(jj(X)  |a>  [2-18al 


a.,(X) 


<g| D.  I fXf\ D^\g>  < g\ | fxf\ D \ g> 


“fi  (©,  + ©) 

fg  ' 


(©,  + ©) 

fg  ' 


[2-18b] 
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where  ajk(x)  is  considered  to  be  an  operator  that  connects  nuclear 


states  1 a> 

and  1 b> . 

When 

to  « Ofg,  which 

is 

the  case  for 

nonresonant 

scattering. 

the 

denominator  may 

be 

treated  as  a 

constant.  The  great  utility  of  Placzek  theory  is  that  it  removes 
the  necessity  of  calculating  excited  electronic  states  of 
molecules . 

These  approximations  permit  attention  to  be  directed  on  the 
motion  of  nuclei,  and  thus  precede  the  harmonic  approximation 
typically  described  in  introductory  texts. 

Spectral  Density  and  Polarizability  Fluctuations 

Before  proceeding  with  the  mathematical  treatment  of  nuclear 
motion,  it  is  useful  to  clarify  its  connection  with  measurements 
made  experimentally.  The  time-correlation  function  approach  that 
follows  also  illustrates  that,  while  Rayleigh  scattering  is  an 
elastic  process  and  Raman  scattering  an  inelastic  one,  they  share 
a common  physical  origin  in  their  dependence  upon  polarizability 
fluctuations  that  arise  from  the  motion  of  nuclei^^. 

Experimental  Geometry  and  General  Equations. 

With  reference  to  Figure  2-1,  consider  a sample  located  at  the 
origin.  Provided  the  scattering  medium  is  nonmagnetic  and 
nonconducting,  the  description  of  its  response  to  an  incident 
electric  field  may  be  given  in  terms  of  a local  polarizability. 
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a(r,t)  = ocl  + boi(r,  t) 


[2-19] 


in  which  I is  the  second  rank  unit  tensor,  and  6a(r,t)  is  the 

polarizability  fluctuation  tensor  at  position  r and  time  t.  By 
convention,  the  polarization  of  the  laser  is  taken  to  be  along  the 
z-axis.  Let  the  incident  electric  field  be  a plane  wave  of  the 
form 


in  which  n^  is  a unit  vector  in  the  direction  of  the  electric 

field,  Eq  is  the  field  amplitude,  is  its  wave  vector  and  ca^  its 
angular  frequency. 

For  an  experimental  arrangment  that  includes  a polarizer 
admitting  scattered  light  with  polarization  n^,  and  a monochrometer 

selecting  light  of  angular  frequency  cOf,  the  spectral  density  of 

light  scattered  into  the  detector  will  be  proportional  to. 


[2-20] 


[2-21] 


CO  = CO 


laser 


CO 


scat 


where  the  brackets  indicate  the  ensemble  average,  and  the 
subscripts  i and  f indicate  the  relevant  initial  and  final 
coordinate  axes  in  the  laboratory  frame.  The  proportionality 
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constant  is  the  product  of  two  factors,  the  first  of  which  depends 
on  the  distance  from  the  sample  to  the  detector,  R,  the  incident 
light  intensity,  Iq  = IEqI^,  and  the  scattering  wave  vector,  k^,  and 
is  given  by. 


0 


16n^R2g,2 


[2-22] 


The  second  factor  is  not  so  easily  obtained  since  it  depends  on 
experimental  details  such  as  the  efficiency  of  the  polarizer  and 
detector. 

It  is  convenient,  for  reasons  which  will  be  made  apparent 
shortly,  to  define  two  quantities  describing  the  intensity  of 
scattered  light  at  the  detector: 


I. 

ISO 


(“)  = -Twl”)  - 


I . (M)  = 

am  SO  ' ' VH 


[2-23] 


which  are  known  as  the  isotropic  and  anisotropic  (or  depolarized) 
spectra.  The  first  subscript  indicates  the  polarization  of  the 
incident  light,  traditionally  labeled  "vertical",  and  the  second 
subscript  indicates  the  polarization  of  scattered  light  at  the 
detector,  which  may  be  arranged  by  the  experimentalist  to  be  either 


vertical  or  horizontal. 
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This  provides  the  general  approach  that  is  appropriate  for 
the  study  of  Raman  scattering  from  a wide  variety  of  molecules,  and 
the  more  difficult  case  of  molecular  liquids. 

Special  Cases  of  Interest 

The  description  of  light  scattering  from  molecular  liquids 
quickly  becomes  complicated  by  molecular  anisotropy  and 
intermolecular  interactions.  Thus  it  is  useful  to  have  in  mind  the 
simpler  results  of  certain  special  cases. 

In  the  simplest  case,  the  medium  consists  of  independent, 
isotropic,  and  structureless  scattering  centers^^.  The 

polarizability  tensor  reduces  to  a scalar,  a = a 5^^,  and  the 

polarizability  fluctuations  assume  a particularly  simple  form: 

5oi.^(q,t)  = (ii/  exp[i  (qr  • (t)  ) ] [2-24] 

where  the  sum  is  performed  over  the  scatterers  within  the  relevant 
volume.  As  a result,  the  (Rayleigh)  scattering  is  determined  only 
by  translational  motion  which,  for  a system  of  identical 
scatterers,  can  be  expressed  in  terms  of  density  fluctuations  in 
the  scattering  volume  (as  indicated  by  the  subscript  "SV") , 
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J^(q-,  co)  = I^(q,oi)  =0 

I^{q,a)  = dt  <5p*^(<i,  0)  5p^^,(<j,  t)  > 

J — OO 

bp^^(q,t)  =j  d^rbp(r,t)  exp[ig-r] 


[2-25a] 

[2-25b] 

[2-25c] 


It  can  be  seen  that  the  anisotropic  scattering  vanishes 
identically,  as  expected.  To  proceed,  one  typically  uses  kinetic 
molecular  theory  or  a hydrodynamic  model  to  calculate  these  density 
fluctuations'*^ . 

The  significance  of  the  depolarized  scattering  becomes 
apparent  if  we  relax  the  requirement  of  molecular  isotropy,  though 
we  will  retain  the  assumption  that  these  molecules  do  not  interact 
with  one  another.  For  such  a system  of  N identical  molecules,  and 
assuming  the  statistical  independence  of  translation  and  rotation, 
the  time  correlation  function  of  interest  is'*^ 

= <N><oi.^(0)*  a.^{t)>F^{q,  t)  [2-26a] 

F (q,t)  =<exp[iq-  (r  (t)  -r  (0))]>  [2-26b] 

* J J 

The  aif(t)  are  elements  of  the  molecular  polarizability  tensor  in 

the  laboratory  fixed  coordinate  system  and  change  with  time  because 
of  molecular  reorientation.  The  q dependence  can  be  seen  to  lie 
entirely  in  density  fluctuations,  formally  separated  as  F^Cq, t). 
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The  tensorial  nature  of  the  polarizability  can  no  longer  be 
avoided,  but  it  can  be  simplified  by  first  choosing  a molecule- 
fixed  coordinate  system  such  that  the  polarizability  assumes 
diagonal  form  and  which  contains,  as  its  components,  the  three 

principal  values  of  the  polarizability,  a^,  tt2,  ttg.  One  may  then 
write  the  tensor  as  the  sum  of  an  isotropic  part  Oq  and  anisotropic 
part  P: 


a = + p 


a 


0 


Tr  [a] 


Tr[p]  = 0 


[2-27a] 

[2-27b,  2-27c] 


In  the  simplest  non-trivial  case,  that  of  a cylindrically 
symmetric  molecule,  one  may  label  the  axis  of  highest  symmetry  by 
II  , and  the  remaining  two  by  ±.  This  gives. 


a 


0 


+ 2a  ) 

± 


[2-28a,  2-28b] 


The  position  of  the  detector  with  respect  to  the  sample  is 
ordinarily  fixed  such  that  the  scattering  vector  q is  fixed  at  90°. 
The  dependence  on  q may  be  then  left  implicit  and  the  spectral 
densities  assume  particularly  simple  forms: 


44 


-Ti^^(co)  = J^(co)  - = <N>a^(co) 


[2-29a] 


I^(co)  =<N>^^^{oi) 


[2-29b] 


Although  this  discussion  has  assumed  cylindrical  symmetry, 
equations  2-29a  and  2-29b  are  general  in  the  sense  that  only  the 
depolarized  spectrum  will  contain  information  about  the  optical 
anisotropy,  however  its  precise  form  will  depend  on  the  arrangement 
of  atoms  within  the  particular  molecule. 

In  the  final  special  case,  the  effects  of  intermolecular 
interaction  will  be  considered.  At  moderate  densities,  even  the 
ideal  gases  argon  and  krypton  have  been  found  to  depolarize 
scattered  light.  The  anisotropy  responsible  for  this 
depolarization  is  a result  of  the  electron  cloud  deformation  that 
occurs  when  two  atoms  approach  one  another.  If  the  collisional 
pair  is  regarded  as  a quasilinear  molecule  for  the  duration  of  the 
collision,  we  may  write  the  polarizability  tensor  in  the  form 


where  e^j  is  a unit  vector  in  the  direction  of  - r j . In  the 

limit  of  infinite  separation,  and  for  identical  atoms,  these  terms 


[2-30] 


reduce  to: 
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liiti 


[2-31a] 


lim  P (r . ,)  = 0 


[2-31b] 


Both  terms  will  reflect  the  effects  of  the  collision.  The 


depolarized  spectrum  will  be  proportional  to  the  Fourier  transform 
of  the  time  correlation  function 


in  which  the  brackets  indicate  the  ensemble  average  and  P2  is  the 
Legendre  polynomial,  P2(x)  = kj  (3x2  _ i),  jf  one  can  neglect  all 
but  two  body  collisions,  then  the  symmetry  of  the  colliding  system 
simplifies  the  result  to 


Equations  2-32  and  2-33  both  omit  a density-dependent 
multiplicative  constant.  These  results  are  the  simplest  examples 
of  what  is  denoted  "Collision  Induced  Light  Scattering",  or  CILS^^. 

The  description  of  light  scattering  from  liquids  containing 
anisotropic  molecules  will  contain  features  in  common  with  each  of 
the  special  cases  described  above  as  well  as  specific  models  for 


[2-33] 
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the  molecular  motion  involved.  These  models  are  conveniently 
classified  according  to  the  modes  of  motion  treated. 

Models  of  Polarizability  Fluctuations 

Vibrational  Motion:  The  Harmonic  Approximation^^ 

If  we  exclude  resonant  light  scattering,  the  type  of  motion 
which  occurs  at  the  fastest  time  scales  and  with  the  largest 
changes  in  energy  is  vibrational  motion.  Vibrational  motion,  which 
must  be  treated  quantum  mechanically,  is  exhaustively  treated  in 
the  literature;  consequently  only  the  most  important  features  will 
be  discussed  here. 

Let  us  assume  that  the  electronic  problem  for  a particular 
molecule  has  been  solved  and  its  equilibrium  geometry  found.  The 
Hamiltonian  describing  vibrational  motion  for  a molecule  consisting 
of  n atoms  oscillating  about  their  equilibrium  positions  can  be 
written  as 


osc 


+ ()) 


[2-34] 


where  Pi^/2mi  represents  the  kinetic  energy  of  the  atom  i and  (j)  is 

the  atomic  potential  which  is  a function  of  the  instantaneous 
positions  of  the  atoms. 
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Provided  the  oscillation  amplitudes  of  the  atoms  are  small,  it 

is  a good  approximation  to  expand  (|)  in  a power  series  of  the 

coordinates  describing  the  deviation  from  the  equilibrium  positions 
as 


[2-35] 


where  the  repeated  indices  over  a and  P indicate  summation  over 

those  indices,  (j>°  is  the  potential  energy  of  the  atoms  at  the 

equilibrium  positions,  and  fab(i,j)  the  generalized  force  constant 
given  by 


d^(\) 
dx.  dx 

la 


a,  P = x,y,  z 


[2-36] 


wwhich  forms  a 3n  x 3n  matrix.  The  quantity  ^{i)  represents  the 

a component  of  the  coordinate  that  represents  the  deviation  of  the 

i’^*'  atom  from  the  equilibrium  position  and  is  known  as  its 
vibrational  coordinate.  The  linear  terms  in  the  vibrational 
coordinates  vanish  because  the  atoms  are  in  their  equilibrium 
positions . 

Since  the  function  (j)°  does  not  depend  on  the  vibrational 


coordinates  it  is  convenient  to  define  a new  Hamiltonian  given  by 
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H 


osc 


OSC 


[2-37] 


This  is  simplified  further  by  introducing  the  mass-reduced 
coordinates  Ua(i)  and  mass-reduced  generalized  force  constants 
9ab(if  j)  which  are  given  by 


u (i)  = ^ (i)m^^^ 

a ' ' 1 


[2-38a] 


(i,  j) 


(m.m.) 

' 1 j' 


-1/2 


[2-38b] 


Expressed  using  these  mass-reduced  variables,  is  given  by. 


H 


OSC 


1 

2 


5r„p(i,  j)  uji)  Up(j)  + 


V 


pji)  - uji)  pf  = E„p„Mi) 


[2-39a] 

[2-39b,  2-39c] 


where  represents  the  enharmonic  term,  i.e.,  the  sum  of  all  those 
higher-order  terms  not  written  out  explicitly;  the  harmonic 
approximation  consists  in  neglecting  this  term. 

One  proceeds  by  finding  a linear  transformation  from  {Ua(i)} 
to  a new  set  of  coordinates  {gj}  such  that  is  expressed  as  a 
linear  combination  of  only  squares  of  the  new  coordinates,  the 
cross-product  terms  having  been  eliminated.  This  is  called  a 
canonical  form  and  the  coordinates  are  known  as  the  normal 
coordinates.  In  these  coordinates  the  Hamiltonian  is  given  by. 
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[2-40] 


This  leads  to  the  well  known  result  that  the  eigenenergy  for  the 
eigenstate  |n,j>  is  given  by  = (h/27t)0|j(n  + h.)  from  which  the 

vibrational  frequency  v = &/2%  is  readily  obtained. 

For  an  ensemble  of  molecules,  the  spectral  density  for  the 
scattered  light  associated  with  transition  of  the  nuclear  state  | a> 
to  |b>,  observed  in  the  j polarization  and  with  the  incident  light 
polarized  in  the  k direction  is  obtained  from  equation  2-13,  which 
is  repeated  here,  in  a slightly  altered  form,  for  convenience: 


4c^ 


< > 5 (CO  - - CO,,) 


0 


ab 


[2-41] 


Here  the  a,b  superscripts  have  been  suppressed  for  brevity,  and  the 
brackets  now  indicate  the  ensemble  average. 

By  introducing  the  Fourier  representation  of  the  delta 
function  equation  2-41  can  be  manipulated  into  a form  explicitly 
indicating  the  relevant  time  correlation  function^®. 


dt  exp[-iAcot]<a^^(0)  *oi.^(t)  > 


Aco  = CO  “ Qq 


C = Ip(coV4c') 


[2-42] 
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It  may  be  observed  that  this  is  formally  equivalent  to  the  result 
given  in  equation  2-21.  The  isotropic  polarizability  time- 
correlation  function  is  given  by 

= -|rr[a]  [2-43] 


For  scattering  from  a fundamental  vibration,  the  isotropic 

I 

Raman  tensor  can  be  expanded  in  terms  of  the  relevant  normal 
coordinate  q giving. 


a 


0 


G. 

ISO 


]ol"<q(0)q(t)> 


[2-44] 


These  are  the  appropriate  formulae  for  light  scattering  from 
a dilute  gas,  but  the  situation  becomes  more  complicated  in 
condensed  phases.  In  general,  the  effects  of  intermolecular 
interaction  on  the  Raman  scattering  spectrum  of  the  normal  mode  can 
be  described  by  the  interaction  Hamiltonian  V given  by 


y = V 


0 


It 


(V.4>q. 

' IJ  ^1 


+ 


J-J 


[2-45] 


where  V°  does  not  depend  on  the  normal  coordinate  and  q^  is  the 
normal  coordinate  associated  with  molecule  i.  The  interaction 
potential  given  here  is  a function  of  normal  coordinates  because  of 
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molecular  vibration  that  modulates  the  intermolecular  potential. 
It  is  of  the  form  applicable  to  the  intermolecular  dipole-dipole 
interaction,  quadrupole-quadrupole  interaction,  hydrogen-bonding 
interaction,  as  well  as  other  types  of  interactions. 

Low  Frequency  Rotational  and  Translational  Motion 

Although  the  quantum  mechanical  treatment  of  rotational  motion 
is  well  known,  a discussion  of  its  results  is  not  appropriate  here, 
as  a glance  at  the  experimental  spectra  will  confirm:  in  molecular 
liquids,  the  inhomogenous  line  broadening  of  the  vibrational  and 
rotational  spectra  by  intermolecular  interactions  makes  the  fine 
rotational  structure  of  the  intramolecular  vibrational  peaks 
indistinct . 

The  strong  intermolecular  interactions  present  in  a liquid 
also  complicate  the  low  frequency  shift  region  so  that  two 
definitions  are  in  order  before  proceeding.  Rayleigh  scattering 
refers  to  the  elastic  and  quasi-elastic  scattering  of  light  by 
density  fluctuations  to  give  a broadened,  but  unshifted,  line 
whereas  Raman  scattering  refers  to  the  inelastic  energy  transfer 
from  the  incident  light  into  (or  from)  a quantized  energy  state  to 
create  a new  photon  frequency  shifted  from  the  incident  light.  In 
dilute  gases  these  regions  are  separated  by  at  least  1000  cm"^  and 
they  are  studied  using  distinct  theoretical  and  experimental 


methods . 
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Under  the  relatively  high-resolution,  low-gain  experimental 
conditions  of  Rayleigh  scattering,  the  Rayleigh  wing  exhibits 
distinct  non-Lorentzian  behavior  which  under  the  low-resolution, 
high-gain  conditions  of  the  Raman  experiment  appears  as  a weak 
local  maximum  on  the  wing  of  the  intense  Rayleigh  line.  These 
"peaks"  have  been  observed  in  essentially  all  isotropic  or 
anisotropic  liquids  studied  to  date.  Though  there  are  important 
differences  in  the  details  of  their  interpretation,  their  is  strong 
evidence  to  suggest  that  for  anisotropic  liquids  the  additional 
depolarized  scattering  has  a predominantly  librational  origin, 
while  the  isotropic  spectrum  contains  effects  from  both  librational 
and  translational  motion.  The  latter  is  often  interpreted  in  terms 
of  intermolecular  vibrational  motion  while  the  former  may  be 
thought  of  as  a quantized  molecular  libration  much  like  those  that 
give  rise  to  rotatory  Raman  lattice  modes  in  the  solid  phase,  or  as 
an  incomplete  loss  of  memory  after  collision  in  the  rotational- 
diffusional  behavior  molecules.  For  spherically  symmetric 
molecules  such  as  CCI4  and  SFg,  Bucaro  and  Litovitz  and  Posch  and 
Litovitz  have  been  able  to  account  for  the  depolarized  spectra  by 
considering  the  polarizability  anisotropy  caused  by  both  collision- 
induced  frame  distortion  and  dipole-induced  dipole  mechanisms. 
Similar  approaches  have  been  successful  in  the  description  of  the 
scattering  from  anisotropic  molecules,  specifically  water. 
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Consequently,  this  discussion  will  cover  a treatment  which  is 
appropriate  for  molecular  liquids  and  whose  formalism  can  be 
readily  extended  to  light  scattering  from  both  the  lower  frequency 
rotational  motion  of  molecules  in  their  ground  vibrational  state 
and  to  their  translational  motion.  The  treatment  relies  heavily  on 
the  insight  provided  by  the  study  of  low  frequency  collision 
induced  light  scattering  from  rare  gasses  which  occurs  at  moderate 
densities  and  is  commonly  denoted  the  Dipole-Induced-Dipole  (DID) 
model'*®. 

The  formalism  begins  by  adopting  a physically  broader 
perspective  which  will,  unfortunately,  require  some  changes  in 
notation.  Consider  an  element  of  volume  containing  a number  of 
scattering  centers.  From  the  preceeding  discussion  we  can 
immediately  write  the  spectral  density  of  the  scattered  light  in 
terms  of  the  polarizability  of  the  scattering  centers, 

HV,n,m)  = E<.,b  [2-46aJ 

(0)>1  [2-46b] 

where  n and  m are  unit  vectors  in  the  direction  of  the  incident  and 
observed  polarization  respectively.  Greek  subscripts  refer  to 
Cartesian  components  in  the  laboratory  frame  and  FT  refers  to  the 
Fourier  transform,  which  is  defined  by 


FT[f(t)]  = -^f°°dt  exp[-2nivt]  f(t) 

2n 
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[2-46c] 


is  the  effective  polarizability  of  the  i*^^  molecule  and  is 


defined  by 


[2-47] 


where  is  the  total  dipole  moment  on  the  molecule,  and  E°  denotes 

the  incident  laser  field.  An  expression  for  the  total  dipole 
moment  can  be  obtained  from  long-range  perturbation  theory^^: 


in  which  refers  to  the  permanent  dipole  moment  on  the  ith 

molecule.  The  remaining  terms  are  refered  to  as  the  "collision 
induced"  contributions  to  the  total  dipole  moment.  They  contain 

contributions  from  the  permanent  molecular  polarizability,  a'^-,  the 
first  hyperpolarizability,  P,  and  the  dipole-quadrupole 


[2-48] 


polarizability,  due  to  the  influence  of  (E)^,  the  local  electric 
field  and  its  gradient,  (E')^.  The  local  field  is  the  effective 
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field  experienced  by  the  i*"^  molecule.  It  is  given  by  the 
superposition  of  the  external  field  with  those  generated  by  the 
multipoles  of  nearby  molecules: 


Eb  *Ei  Ebv  tc’lijjei  • ... 


[2-49] 


is  the  n-pole  propagator  defined  by. 


T (n) 


(ij)  = (-i)"[v^  v„ . . .y 


“l  “2 


Irl  ^ 


[2-50] 


0^  is  the  effective  quadrupole  moment  on  the  ith  molecule,  defined 
by. 


= 

ap 


[2-51] 


The  final  expression  for  the  effective  polarizability,  is 
given  in  Cartesian  components  to  first  order  in  the  n-pole 
propagator  and  written  assuming  the  Cartesian  summation 
convention,  i.e.,  summation  is  implied  over  repeated  subscripts: 


P ^ 

aP 


— 


= a 


«3 


n\  = 

ap 


a 


aP 


y^.  . [«■“ 


— T ^ 

2 ore 


( ij  ) Vy5  3 Vy5  ^ ^ J ) «eP+  * ' ' ^ 


[2-52] 


It  has  been  written  to  indicate  the  separate  contributions  of  the 
permanent  molecular  polarizability,  and  an  induced  part,  11^.  The 
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resulting  expressions  for  the  isotropic  and  depolarized  spectral 
densities  are, 

-I 

J.  (v)  = < \FT  ( V.  —52  P"  (t)  ) > 

[2-53] 

J . (v)  = <\FT['V.  —Y'  . [Piit) 

am  SO  • ^ Q it— /a3  a3 

, [2-54] 

- 5 „ (—  52  PMt)  ) ] ] p > 

a3  2 YY  ' J J I 

Chapter  4 contains  a more  detailed  description  of  how  these 
terms  may  be  calculated  from  configurations  generated  by  a 
molecular  dynamics  simulation.  The  procedure  requires  using  the 
well-known  relationship  between  properties  given  in  the  laboratory 
coordinate  system,  labeled  by  Greek  subscripts,  and  coordinates 
fixed  in  the  molecular  frame,  labeled  by  Latin  subscripts.  Thus, 
one  may  calculate  the  contribution  due  to  the  permanent  molecular 
polarizability  using 

= Ei  Ej  C i Cp,  [2-55] 

where  C is  the  linear  transform  specifying  the  orientation  of  the 
two  coordinate  systems  with  respect  to  each  other. 

Some  general  comments  on  the  interpretation  of  these  spectral 
densities  are  in  order  here  before  postponing  a discussion  of  the 
remaining  details  until  Chapter  4 . The  time  correlation  function 
of  the  total  polarizability  tensor  can  be  separated  into  a sum  of 


57 


terms  containing  only  self-correlations  of  the  permanent  molecular 
polarizabilities  and  a sum  over  terms  containg  cross-correlations. 
The  so-called  "allowed  spectrum"  is  the  Fourier  transform  of  the 
sum  over  self-correlations.  It  is  related  to  single  particle 
translational  motion  through  the  trace  of  the  permanent  molecular 

polarizability,  oLq,  and  to  orientational  motion  through  the 

anisotropy,  5,  of  the  permanent  molecular  polarizability. 

The  remaining  terms  constitute  the  Collision  induced  spectrum. 
The  "PTH"  and  the  "aTA"  are  the  most  important  of  these  terms.  The 

former  includes  terms  involving  "aTa",  which  is  due  to  isotropic 
part  of  polarizability  alone  and  is  relatively  large,  and  two 
lesser  contributions,  "aT6"  and  "6t6",  which  are  much  smaller  and 
are  due,  in  part,  to  the  optical  anisotropy. 

Reduction  Procedures  and  Relationships  to  Other  Methods 

The  dynamical  region  probed  by  spontaneous  light  scattering 
overlaps  those  probed  by  other  techniques,  particularly  infrared 
(IR)  spectroscopy,  and  inelastic  neutron  scattering  (INS) . This  is 
one  of  the  motives  for  recasting  the  expressions  for  the  scattered 
intensity  into  a form  suitable  for  comparison  with  these  other 
techniques.  The  second  motive  is  the  desire  to  distinguish  the  low 
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frequency  inelastic  processes  from  the  much  more  intense 
contribution  of  elastic  or  Rayleigh  scattering'*®. 

The  Raman  spectral  density,  I (©) , is  given  by  the  Fourier 

transform  of  the  total  polarizability  time  correlation  function  and 
can  be  expressed  as  the  sum  of  the  intensities  on  the  Stokes  side, 

(©)  , and  the  anti-Stokes  side,  I®®(©): 


where  m is  the  shift  in  angular  frequency.  A function  R((B)  can  be 
defined  according  to. 


J(cd)  = + J^®(co) 


[2-56] 


R(co)  = co[  + J^®(co)  ] 


[2-57] 


According  to  the  principle  of  detailed  balance 


= exp[-hQ/2nkr] 


[2-58] 


which,  on  substitution  into  equation  2-57,  gives  us: 


R(co)  = 0l){  1 -exp  [ -hco/2nkT]  } J®‘^(co) 


[2-59] 


Taking  into  account  the  dependence  of  the  scattered  intensity  on 
the  fourth  power  of  the  real  frequency  the  final  expression  for 
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what  is  commonly  termed  the  "reduced"  representation  of  the 
spectrum  is. 


where  Qq  is  the  frequency  of  the  laser.  Typically,  the  angular 

frequency  shift  is  converted  to  the  wavenumber  shift.  The  R(©) 

function  has  been  interpreted  as  the  energy  absorbed  in  the 
scattering  process. 

In  one  common  variation,  the  reduced  representation  is 
reported  in  the  mathematically  equivalent  form: 


In  the  low-frequency  region  the  relative  intensity  dependence 
upon  the  fourth-power  term  is  not  important.  For  example,  using  an 
Argon-ion  laser  at  514.5  nm  as  the  exciting  source,  as  is  the  case 
for  the  experimental  results  reported  in  this  work,  the  influence 
of  this  factor  in  the  region  10  to  200  cm‘^  is  around  4%. 
Neglecting  the  fourth-power  term  and  expanding  the  exponential, 

keeping  on  the  first  term  linear  in  to,  gives  an  approximate 
expression  for  R(©)  that  is  in  common  use: 


R(co)  = (cOg  - co)  ^ co{l  - exp  [ -hco/2nkr]  } 


[2-60] 


R(q)  = J®'^(co)  (cOq  - co)  “*co  tanh  (hco/4nkT) 


[2-61] 


approx  (^)  „ ^2 


[2-62] 
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Ordinarily,  one  measures  only  the  Stokes  intensity  and  the 
superscript  is  suppressed. 

The  R(©)  representation  is  convenient  because  it  permits 

direct  comparison  with  the  infrared  absorption  coefficient,  aij.(co), 
which  has  a similar  relationship  to  the  dipole  correlation 
function,  <|0,  (0)  |j,  (t)  >: 

a.^(a)  « “ tanh  (hco/4nkT)  [2-63a] 

= /*“  dt  <p(0)p(t)>  [2-63b] 

J “OO 

The  reduced  representation  is  also  useful  for  coping  with  the 
enormous  intensity  of  the  central  Rayleigh  line  that  is  a serious 
problem  in  studies  of  low-frequency  Raman  spectra.  The  time  scale 
for  motions  giving  rise  to  the  Rayleigh  line  is  relatively  slow 
compared  to  events  reflected  in  the  low-frequency  Raman  region. 
The  correlation  function  can  be  expressed  as  an  exponential  decay 
giving  rise  to  a Lorentzian  shape  for  the  Rayleigh  line.  A 
Lorentzian  band  centered  at  a Raman  shift  of  zero  cm"^  is  converted 

into  a rather  smooth  baseline  in  the  R(v)  representation,  thus 
making  the  weaker  Raman  features  much  more  apparent. 


CHAPTER  3 

EXPERIMENTAL  DESIGN  AND  RESULTS 

Raman  Cell 

A cell  appropriate  for  light  scattering  at  elevated 
pressures  was  specially  designed  and  constructed.  The  cell 
consisted  of  a series  of  concentric  cylindrical  sleeves  that 
were  held  in  place  by  chamber  heads  at  each  end.  Window 
flanges  at  the  two  ends  completed  the  basic  assembly.  The 
innermost  sleeve  was  constructed  of  quartz  and  had  an  inner 
diameter  of  0.240  inches,  an  outer  diameter  of  0.475  inches, 
and  a length  of  2.495  inches;  the  sample  was  contained  by  this 
sleeve.  A brass  sleeve  (0.484  inches  i.d.,  0.625  inches  o.d.) 
surrounded  the  first  sleeve  and  had  an  overall  length  of  2.031 
inches.  A shoulder  extended  radially  from  0.610  inches  from 
center  and  was  0.065  inches  deep  from  either  end.  Three  slots 
were  cut  in  this  sleeve,  each  1/8  inch  wide  and  1.125  inches 
long  centered  longitudinally  with  their  centers  90  degrees 
apart . 

The  brass  sleeve  was  followed  by  a second  quartz  tube 
(0.870  inches  i.d.,  0.988  inches  o.d.,  2.093  inches  in 
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length) . Coolant  fluid  (50%  ethylene  glycol  + 50%  water  by 
mass)  flowed  through  the  gap  between  the  brass  and  outer 
quartz  sleeve.  A final  sleeve  was  constructed  of  delrin 
(1.103  inches  i.d.,  1.312  inches  o.d.)  and  protected  the 
bundle  from  outside  shock.  It  had  an  overall  length  of  1.956 
inches  and  a shoulder  extending  0.306  inches  from  either  end. 
Three  slots  were  cut  into  the  delrin  sleeve  to  line  up  with 
those  in  the  brass  sleeve.  The  slots  were  3/8  inches  wide, 
and  1.3125  inches  long,  again  centered  longitudinally. 

The  chamber  heads  were  designed  to  permit  high  pressures 
within  the  sample  container  while  avoiding  compression  stress 
that  could  shatter  the  quartz  components.  This  was 
accomplished  by  machining  a series  of  concentric  shoulders  on 
the  side  facing  the  sleeves;  they  permitted  the  quartz  sleeves 
to  float  between  the  two  chamber  heads  without  breaking  the 
seals.  The  chamber  heads  were  constructed  of  stainless  steel 
and  machined  to  a thickness  of  0.75  inches.  The  central  hole 
was  drilled  through  to  a diameter  of  0.245  inches.  The 
innermost  shoulder  extended  radially  from  the  central  hole  to 
a diameter  of  0.484  inches  and  was  3/8  inches  deep.  There 
followed  a shoulder  extending  radially  from  the  first  to  0.614 
inches  from  the  center  of  the  piece.  Finally  a groove  for  the 
outer  quartz  and  delrin  sleeves  was  cut  to  a depth  of  0.120 
inches;  this  groove  had  an  inner  diameter  of  0.854  inches  and 


outer  diameter  of  1.106  inches. 
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Figure  3-1.  Chamber  head  component  of  Raman  Cell. 
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Figure  3-2.  Window  flange  component  of  Raman  scattering  cell. 
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Figure  3-4.  Quartz  sleeves. 
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Figure  3-5.  Cut-away  view  of  assembly. 
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Figure  3-8.  Block  diagram  of  experimental  setup. 
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Figure  3-9.  Calibration  of  pressure  gage  using  methane  gas. 
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The  opposite  (outside)  face  of  the  chamber  head  had  a 5/8 
inch  central  hole  to  house  the  quartz  window.  0-rings  placed 
in  grooves  at  the  base  of  the  window  housing  sealed  the  window 
against  the  chamber  head. 

Each  chamber  head  contains  two  inlets,  one  for 
introduction  of  the  sample  and  the  other  for  circulation  of 
coolant.  The  stainless  steel  tubes  (1/8  inch  o.d.  for  the 
sample  inlet  and  1/4  inch  o.d.  for  coolant)  were  brazed  into 
the  chamber  head.  The  window  flange  was  held  to  the  chamber 
head  by  six  6-32  bolts  through  holes  drilled  on  a 0.750  inch 
diameter  bolt  circle. 

The  components  of  the  Raman  cell  are  illustrated  in 
Figures  3-1  through  3-4  while  their  assembly  is  illustrated  in 
Figure  3-5.  In  Figure  3-6  the  heat  exchanger  is  illustrated 
which  cooled  gaseous  samples  before  their  introduction  into 
the  Raman  cell  and  Figure  3-7  illustrates  the  one-way  pump 
used  to  stir  the  sample. 

Instrumentation  and  Computer  Interface 

A block  diagram  of  the  assembled  components  is  given  in 
Figure  3-8.  Source  radiation  was  provided  by  an  Argon  ion 
laser  (Spectra  Physics  Model  2020)  tuned  to  the  514.5  nm  line. 
The  power  of  the  incident  laser  beam  was  measured  using  a 
Scientech,  Inc.,  Model  361  powermeter  at  200-400  mW.  The 
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excitation  power  employed  was  300  mW  and  varied  by  no  more 
than  5%.  After  the  initial  polarization  was  confirmed,  only 
one  polarizer  was  used  and  it  was  placed  between  the  sample 
and  the  entrance  to  the  monochrometer. 

The  monochrometer  employed  was  a Spex  model  1870  driven 
by  a Spex  model  1673c  minidrive  at  the  rate  of  0.05  nm/sec 
which  permitted  acquisition  of  the  4000  wavenumber  spectra  in 

approximately  8 hours.  A slit  width  of  150  |xm  was  used 

providing  a spectral  resolution  of  approximately  4 cm”’^.  Two 
standards  were  used  to  calibrate  the  monochrometer:  the 
mercury  line  from  fluorescent  room  lights  (1123.4  cm"^)  and  the 
sulfate  peak  at  980  cm"^  of  a 0.50  M S04^“  solution  prepared 
from  K2SO4. 

The  signal  was  detected  with  a single-channel  photon- 
counting system  (RCA  31034A  PMT)  connected  to  a Keithley  195A 
digital  voltmeter.  The  output  was  then  passed  to  a Gateway 
2000  microcomputer  containing  a 486  CPU  operating  at  66 
megahertz  with  32  megabytes  of  RAM.  The  interfacing  protocol 
used  for  the  connection  of  the  electronic  hardware  to  the 
micro-computer  was  the  standardized  IEEE  general  purpose  bus 
interface  (GPIB) , including  the  8-bit  i/o  card  and  its 
accompanying  software. 

Pressure  data  were  acquired  with  a Sensotec  model 
TJE/7313-22  strain  gage  pressure  transducer  which  has  an 
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effective  pressure  range  of  0-1,000  psi.  This  pressure 
transducer  was  used  together  with  a Beckman  Industrial  600 
Series  electronic  readout  device.  The  Sensotec  transducer  was 
calibrated  against  a Ruska  Model  2465  standardized  Dead  Weight 
Pressure  Gage.  The  Ruska  gage  had  a pressure  range  of  0.000 
psia  to  650.000  psia  with  an  accuracy  of  ±0.001  psia.  The 
pressure  calibration  results  are  shown  in  Figure  3-9. 

A two-wire  platinum  resistance  temperature  device  (RTD) 
was  used  in  the  experiments.  The  ceramic,  cylindrically- 
shaped  RTD  was  used  in  conjunction  with  a Keithley  195A 
digital  multimeter  to  read  resistance  data.  The  resistance 
values  were  converted  to  their  corresponding  temperatures,  in 
Celsius,  using  the  following  equations  and  parameters  in  which 

a = 0.00385,  5 = 1.502,  Rq  = 100.00  O,  and  R^  is  the 

resistance  of  interest. 
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Figure  3-11.  The  three  normal  vibrational  modes  of  water. 
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Figure  3-13.  Reduced  Spectrum  of  H20(l)  resolved  into  its  components. 
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Figure  3-14.  Low  wavenumber  depolarized  Raman 
scattering  spectrum  of  H20(l)  at  2 °C. 
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Figure  3-15.  Low  wavenumber  depolarized  Raman  spectrum 
of  H20(1)  at  8 °C. 
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Figure  3-16.  Low  wavenumber  depolarized  Raman  spectrum 
of  H20(1)  at  15  °C. 
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Figure  3-17.  Low  wavenumber  depolarized  Raman  spectrum 
of  H20(1)  at  25  °C. 
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The  pressure  and  temperature  measuring  devices  were 
interfaced  to  the  micro-computer  via  the  IEEE  general  purpose 
interface  bus.  All  data  were  acquired  with  computer  programs 
written  specifically  for  these  experiments  in  ANSI  C standard 
code . 

Temperature  control  was  performed  by  a Fisher  Scientific 
Chiller  capable  of  controlling  the  coolants  temperature  to 
±0.1  K.  All  actual  temperature  measurements  were  performed 
using  the  RTD  immersed  in  the  sample. 

Gases  and  Chemicals 

Doubly-distilled  water  from  a MilliPore  Milli-Q  UV  Plus 
UltraPure  Water  System  was  used  after  degassing.  Methane  and 
gas  was  obtained  from  Scott  Specialty  Gases  and  had  a purity 
greater  than  99.99  mole  percent.  The  carbon  dioxide  gas  was 
also  obtained  from  Scott  Specialty  Gases  and  also  had  a purity 
greater  than  99.9  mole  percent.  Both  were  used  without 
further  treatment.  Tetrahydrofuran  (spectroscopic  grade)  was 
obtained  from  Aldrich  Chemical  Company  and  used  after 


degassing . 
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Representative  Spectra 
Light  Scattering  from  Water 

Because  clathrate  hydrates  can  be  regarded  as  a distorted 
form  of  ice,  liquid  water  and  ice  serve  as  particularly 
important  benchmark  systems.  Furthermore,  several 
spectroscopic  studies  of  these  systems  are  present  in  the 
literature  which  enabled  confirmation  that  the  Raman  apparatus 
was  performing  adequately. 

Figure  3-10  contains  the  light  scattering  spectrum  of 
liquid  water  at  25  °C  and  it  is  decomposed  into  its  polarized 
and  depolarized  components  in  Figure  3-12.  Some  general 
observations  may  be  made  though  a more  detailed  discussion 
will  be  reserved  for  Chapter  5. 

The  large  wavenumber  region  of  the  spectrum  may  be 
attributed  primarily  to  intramolecular  vibrations  perturbed  by 
intermolecular  coupling.  The  water  monomer  has  three  normal 
modes  of  vibration,  which  are  illustrated  in  Figure  3-12. 
These  are:  the  H-O-H  bending  (occuring  at  approximately  1595 
cm"^  in  the  gas  phase) , the  assymetric  stretch  in  which  proton 
motions  are  out  of  phase  (3756  cm‘^  in  the  gas  phase)  , and 
symmetric  stretching  motion  in  which  the  protons  move  in  phase 
(3562  cm~^  in  the  gas  phase) . All  three  modes  are  broadened 
considerably  and  shifted  by  intermolecular  interactions  in  the 
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liquid  phase:  The  region  from  1625  - 1644  cm"^  corresponds  to 
the  H-O-H  bending  motion  while  the  large,  broad  peak  in  the 
3000-3800  cm"^  region  corresponds  to  the  two  types  of  OH 
stretching  motion.  When  the  data  for  this  peak  are  plotted 
over  a smaller  wavenumber  range,  this  peak  can  be  seen  to 
contain  structure  indicating  the  strength  and  type  of 
intermolecular  interactions  present  in  the  liquid  and  the 
analysis  of  this  spectral  region  along  these  lines  is  an 
important  part  of  Chapter  5. 

The  wide  wavenumber  range  in  these  figures,  together  with 
the  intense  and  broad  Rayleigh  scattering,  tends  to  obscure 
the  features  in  the  low  wavenumber  region.  Consequently,  this 
section,  decomposed  into  its  components,  is  plotted  as  the 
reduced  intensity  in  Figure  3-13.  The  broad,  nearly 
featureless,  increase  in  intensity  extending  from  around  400 
cm"^  to  1000  cm"^  can  be  attributed  to  the  librational  or 
restricted  rotational  motions  of  the  H2O  molecules. 

Walrafen' s classic  review  of  Raman  scattering  in  water 
demonstrated  that  this  region  can  be  decomposed  into  three 
broad  Gaussian  components  centered  near  425,  550,  and  740  cm"^. 
The  three  components  loosely  correspond  to  the  rotational 
modes  of  the  monomer  though  they  are  strongly  perturbed  by  the 
restraints  imposed  by  hydrogen  bonds  present  in  the  liquid 


phase . 
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Hydrogen  bond  stretching  and  bending  motions  give  rise  to 
additional  weak  bands  below  300  cm‘^.  These  may  be  regarded 
as  restricted  translational  modes  of  motion  since  they 
primarily  involve  movement  of  the  much  heavier  oxygen  atoms. 
There  is  a distinct  peak  centered  near  200  cm'^,  corresponding 
to  the  stretching  motions  of  0-H  units,  and  one  that  is  much 
less  obvious  near  90  cm"^  produced  by  the  bending  of  such 
units . 

Considerable  controversy  still  exists  over  the 
interpretation  of  features  closer  to  the  Rayleigh  line. 
However,  we  may  safely  avoid  this  region  by  noting  that  the 
experimental  results  reported  here  do  not  permit  the  extremely 
high  resolution  distinctions  that  will  be  required  to  resolve 
the  debate . 

The  importance  of  the  low  wavenumber  region  lies 
primarily  in  the  information  it  contains  concerning  the 
restrictions  placed  by  intermolecular  forces  on  rotational  and 
translational  motion.  Representative  spectra  showing  the 
temperature  dependence  of  the  low  wavenumber  depolarized  light 
scattering  are  presented  in  Figures  3-14  through  3-17.  The 
decrease  in  intensity  with  temperature  has  been  interpreted  in 
terms  of  a breakdown  of  hydrogen  bonding. 
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The  Water/  Tetrahvdrofuran  System 

The  other  important  benchmark,  for  a particular  hydrate 
forming  system,  is  the  light  scattering  spectra  of  the  pure 
former.  Figure  3-18  displays  the  spectrum  for  pure 
tetrahydrofuran  (THF)  liquid.  The  high  wavenumber  region 
again  contains  features  attributable  to  intramolecular 
vibrations.  C-O-C  ring  stretch  is  responsible  for  the 
spectral  features  occuring  between  850  and  920  cm"^. 
Specifically,  the  sharp  peak  at  913  cm"^  is  due  to  the  C-O-C 
symmetric  stretch. 

Various  CH2  bend  motions  are  responsible  for  the  complex 
of  peaks  that  appear  between  940  and  1600  cm"^.  The  intense 
peak  at  2800  cm'^  is  the  result  of  CH  stretching  motions. 

A representative  spectrum  of  an  aqueous  solution  of  THF 
is  given  in  Figure  3-19.  At  first  glance  it  appears  very 
much  like  simply  the  superposition  of  the  pure  water  and  pure 
THF  spectra,  which  is  not  suprising  since  THF  is  completely 
miscible  in  water.  However,  there  are  subtle  changes  changes 
in  the  CH2  bending  and  CH  stretch  regions  which  become  more 
apparent  when  they  are  plotted  over  a narrower  waveniamber 
range  and  analyzed  quantitatively.  The  changes  in  the  water 
OH  stretching  regions  are  rather  obvious,  indicating  changes 


in  intermolecular  interaction. 
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Figures  3-19  through  3-21  are  included  to  show  the 
dependence  of  the  spectra  on  solution  composition  and  the 
temperature  dependence  of  these  spectra  are  indicated  by 
Figures  3-21  through  3-23.  These  may  be  compared  to  the 
spectra  of  the  solid  clathrate  presented  in  Figures  3-24  and 
3-25.  Figures  3-26  and  3-27  shows  the  decomposition  of  a pure 
THF  spectrum  according  to  polarization.  The  analogous  spectra 
for  the  THF/water  system  is  given  in  Figures  3-28  and  3-29. 
It  will  be  argued  in  Chapter  5 that  a quantitative  analysis  of 
these  spectra  reveal  a great  deal  of  information  concerning 
the  changes  in  intermolecular  interactions  that  are  the  result 
of  changes  in  the  composition  and  temperature  of  the  system. 
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Figure  3-18.  Scattering  from  pure  THF(l)  at  25 
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Figure  3-19.  Light  scattering  from  THF/H2O  solution  at  25  °C  (0.11  mole  fraction  THF)  . 
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Figure  3-20.  Light  scattering  from  THF/H2O  solution  at  25  °C  (0.056  mole  fraction  THF)  . 
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Figure  3-21.  Light  scattering  from  THF/H2O  solution  at  25  °C  (0.04  mole  fraction  THF)  . 
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Figure  3-22.  Light  scattering  from  THF/H2O  solution  at  10  °C  (0.056  mole  fraction 
THF)  . 
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Figure  3-23.  Light  scattering  from  THF/H2O  solution  at  6 °C  (0.056  mole  fraction  THF) . 
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Figure  3-24.  Light  scattering  from  THF  clathrate  hydrate  at 
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Figure  3-25.  Light  scattering  from  THF  clathrate  hydrate  at 
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Figure  3-26.  Depolarized  spectrum  of  THF  at  25 
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Figure  3-27.  Polarized  spectrum  of  THF  at  25 


800 


97 


o o o o 
o o o o 

CO  uo 


Aq.Tsu0q.u1 


-p 

4-1 

"H 

x: 

U2 

<D 

n 

13 

jq 

<D 

> 

(0 

is 


0) 
I — I 

o 

g 

CD 

LO 

o 

o 


u 

o 

LO 

CM 

-P 

g 

o 

-p 

p 

I — I 

o 

CO 

o 

CM 

Pl-I 

H 

4H 

o 

g 

p 

p 

-p 

u 

Q) 

a 

CO 

TJ 

(U 

N 

•H 

P 

cc3 

O 

a 

(U 

p 


• PLJ 

00  m 

CM  E-h 

1 

CO  p 

o 

CU  *H 
P -P 
P o 
tn  fd 
•H  P 
tj  4-1 


98 


Aq.TSueq.ui 


(U 


O 

e 


KD 

LO 

o 

o 


u 

o 

LO 

CM 

-P 

fd 

p 

o 

-H 

-P 

P 

I — I 

o 

in 

O 

CM 

Ijm 

K 

H 

M-l 

O 


P 

-P 

U 

0) 

a 

CO 

T5 
CD 
tSI 
-H 
P 
fd 
• — I 

o 

a 

(U 

Q 


CNJ  Eh 

I 

cn  p 

O 

0)  -H 
p -P 

P u 
tr»  cd 
•H  P 
Pi-i  M-l 


CHAPTER  4 

MOLECULAR  DYNAMICS  SIMULATIONS 
Software  and  Algorithms 

All  molecular  dynamics  (MD)  simulations  were  performed 
using  modules  from  the  TINKER^®  package  of  source  code.  After 
some  modifications,  the  code  was  compiled  and  run  under  the 
Linux  operating  system  on  a pentium  PC.  Although  some 
features  of  the  software  package  are  incomplete,  those 
features  primarily  pertain  to  the  simulation  of  large 
molecules  of  biological  significance  and  were  not  required  for 
the  studies  reported  here;  the  modularity  of  the  code  permits 
one  to  avoid  those  features  entirely.  A user  interface  was 
written  using  C and  the  Tk/Tcl  toolkit  and  all  analysis 
programs  were  written  by  this  author  in  C. 

The  equations  of  motion  for  the  molecular  centers  of  mass 
and  their  orientations  (quaternions)  were  integrated  using  the 
velocity  Verlet  algorithm^®  The  electrostatic  forces  and 
energies  were  evaluated  by  the  Ewald  method^®  while  Lennard- 
Jones  interactions  were  set  equal  to  zero  for  intermolecular 
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separations  greater  than  half  the  length  of  the  simulation 
cell . 


Intermolecular  Potentials 

A wide  variety  of  intermolecular  potentials  have  been 
proposed  to  model  water,  and  the  choice  of  one  requires  some 
defense.  The  simulations  reported  here  used  the  simple  point 
charge  (SPC)  model  introduced  by  Berendsen^^.  This  is  a semi- 
empirical  model  with  parameters  determined  by  fitting  the 
internal  energy  and  pressure  of  liquid  water  to  the 
experimental  values  at  300  K,  and  with  the  further  restriction 
that  the  model  adequately  approximate  the  measured  oxygen- 
oxygen  distribution. 

The  simplicity  of  the  model  makes  it  attractive  from  a 
computational  standpoint  and  the  SPC  model  is  one  of  only  a 
handful  that  have  been  found  to  adequately  describe  the  static 
properties  of  different  polymorphic  forms  of  ice.^^ 
Furthermore,  the  model  has  been  successfully  used  by  others  to 
simulate  solid  clathrate  hydrates^^. 

The  SPC  model  consists  of  three  point  charges  placed  on 
the  atomic  centers  with  the  oxygen  centers  interacting  through 
a Lennard- Jones  potential;  the  water  molecule  is  rigid  and  no 
account  is  taken  of  any  short-range  H-H  or  0-H  forces. 
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Water-guest  and  guest-guest  interactions  are  treated  in 
a similar  manner  with  each  guest  treated  as  a rigid  molecule 
possessing  a number  of  Lennard- Jones  and  electrostatic 
interaction  sites.  Ethylene  oxide  was  represented  by  three 
interaction  sites;  two  methylene  groups  centered  on  the  carbon 
atoms  and  the  etheric  oxygen  atom.  Electrostatic  parameters 
are  those  reported  by  Tse  et  al^^  and  were  evaluated  from  the 
known  dipole  moment  and  structure  of  the  molecule.  The 
Lennard- Jones  interaction  parameters  were  obtained  from  the 
same  source  and  were  those  used  in  the  TIPS2  model  of 
Jorgensen. 

Methane  was  represented  by  five  coulombic  interaction 
sites  and  a single  Lennard-Jones  site;  the  relevant  parameters 
are  also  from  Tse  et  al  and  were  estimated  from  data  on 
krypton  which  has  a solid  state  lattice  spacing  and  binding 
energy  similar  to  methane.  Guest-water  interaction  parameters 
were  obtained  using  the  Lorenz-Berthelot  mixing  rules 
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The  intermolecular  potential  between  two  molecules  m and 


is  given  by 
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and  all  parameters  are  reported  together  in  Table  4-1. 


Table  4-1.  Intermolecular  Potential  Parameters 

Lennard-Jones  Parameters 


Interaction 

0/  A 

4e  /kJ  mol”^ 

Water  0-0 

3.16 

2.60 

Ethylene  Oxide  0-0 

3.05 

3.27 

Ethylene  Oxide 

3.98 

1.91 

CH2-CH2 

Methane  C-C 

3.64 

5.46 

Methane-Water  CH4-0 

3.33 

3.05 

Electrostatic  Parameters 


Site 

Charge  / e 

Water  0 

1 

0 

• 

00 

N) 

Water  H 

0.41 

Ethylene  Oxide  0 

-0.32 

Ethylene  Oxide  CH2 

0.16 

Methane  C 

-0.56 

Methane  H 

0.14 
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Table  4-2.  Asyinmetric  Residue  Unit  for  Host  Lattice  of  SI 
Clathrate  Hydrates^^. 

Lattice  Parameter:  12.0  A (cubic) 

Space  Group:  Pm3n 


Atom 

X 

y 

z 

Population 

Oi 

0.18375 

0.18375 

0.18375 

1 

0, 

o 

• 

o 

0.30822 

0.1173 

1 

Oc 

o 

• 

o 

0.5000 

0.2500 

1 

Hii 

0.2323 

0.2323 

0.2323 

0.5 

Hck 

o 

• 

o 

0.43404 

0.2015 

0.5 

H,c 

o 

• 

o 

0.37874 

0.1588 

0.5 

o 

• 

o 

0.31614 

0.0357 

0.5 

-H 

0.0676 

0.26483 

0.1397 

0.5 

Hih 

0.1170 

0.22763 

0.1588 

0.5 

Oxygen 
is  the 

atom 

site 

indices  c,  i 
of  a proton 

, k refer 
which  is 

to  the  symmetry  sites. 
covalently  bonded  to  0^  and 

hydrogen-bonded  to  0^. 

Notation  from:  International  Tables  for  X-ray  Crystallography 

(Kynoch  Press,  Birmingham  England,  1965),  Vol.  I. 
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Cluster  Simulation  Details 

Simulations  were  performed  for  water  clusters  of  various 
sizes  near  room  temperature.  The  clusters  were  not  confined 
and  periodic  boundary  conditions  were  not  imposed.  Thus,  they 
were  free  to  dissociate  though  this  did  not  occur  in  any  of 
the  simulations  reported  here. 

Initial  configurations  for  the  smaller  clusters  (N=2,  5, 
8)  were  chosen  from  published  reports  of  local  minima 
configurations.  The  larger  clusters  (N=20,  50)  were  arranged 
as  fragments  of  hexagonal  ice.  It  should  be  noted  that  the 
spectra  calculated  here  have  been  generated  using  a single, 
though  long,  MD  run;  thus  it  is  possible  that  the  lower- 
temperature  clusters  may  not  have  had  enough  energy  to 
isomerize  and  may  not  represent  the  complete  distribution  of 
cluster  geometries. 

The  atomic  configurations  so  obtained  were  then  adjusted 
to  the  HOH  bond  angle  and  OH  bond  length  required  by  the  SPC 
model.  The  simulated  sample  was  then  heated  by  uniformly 
increasing  the  velocities  of  each  molecule  to  produce  a 
temperature  increase  of  0.02  K fs"^  followed  by  a 10,  000  MD 
step  equilibration  period  (10  ps)  every  20  K.  This  sequence 
was  repeated  until  the  desired  temperature  was  obtained 
whereupon  the  system  was  again  allowed  to  equilibrate  for 
10,000  steps.  The  dynamics  were  then  followed  for  32768  steps 


(~30  ps)  by  saving  configurations  every  eight  steps.  Cluster 
temperatures  were  calculated  according  to. 


1 2^ 
— m.v  > 

2 ^ ^ 


[4-4] 


Unlike  the  bulk  simulations,  the  clusters  were  permitted 
to  evolve  without  thermal  interference  and  may  be  considered 
constant  energy  systems,  i.e.,  while  the  initial  temperature 
was  deliberately  chosen,  it  was  not  maintained  during  the 
simulation. 


Bulk  Simulation  Details 

MD  simulations  were  performed  for  Canonical  ensembles 
(constant  NVT)  of  water,  ethylene  oxide/water,  and 
methane/water.  Random  initial  configurations  were  used  for  the 
water  simulations,  while  the  SI  x-ray  results  of  Jeffrey  et 
al^^  were  used  to  obtain  initial  configurations  for  the 
ethyelene  oxide  and  methane  simulations.  Asymmetric  residue 
units  are  reported  in  Table  4-2  and  the  discus^^  software 
package  was  used  to  generate  the  full  unit  cell.  The  excess 
protons  were  then  removed  until  the  resulting  configuration 
satisfied  the  correct  valence  of  oxygen  atoms  and  the  correct 
proton  placement  according  to  the  Bernal-Fowler  ice  rules^^. 
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Each  sample  for  the  bulk  water  simulations  consisted  of 
256  molecules  in  a cubic  simulation  cell,  19.46  A on  edge. 
The  methane/  water  systems  were  simulated  by  cubic  cells  24.6 

o 

A on  edge  corresponding  to  eight  SI  unit  cells;  they  contained 
368  water  molecules  and  64  guest  molecules  which  corresponds 
to  the  ideal  clathrate  stoichiometry.  The  initial 

configurations  for  the  ethylene  oxide/  water  simulations  are 
identical  to  those  used  for  methane  with  the  important 
exception  that  the  guests  are  permitted  to  occupy  only  the 
larger  cages  as  experimental  studies  have  indicated. 

Bulk  simulations  of  the  solid  hydrates  followed  the  same 
procedure,  but  with  constant  temperature  conditions  imposed 
during  the  data  collection  steps  using  the  Gronigen  method. 
These  simulations  were  performed  primarily  for  comparison  with 
simulations  performed  by  others. 

Bulk  simulations  of  the  liquids  began  with  a somewhat 
different  heating  procedure.  The  initial  configurations  were 
first  heated  by  the  same  method  used  for  the  clusters  to  400 
K and  equilibrated  there  for  10,000  steps  to  ensure  melting. 
They  were  then  cooled  to  the  desired  temperature  and 
equilibrated  further  for  10,000  steps  before  the  32768  data 
collecting  steps  were  performed  during  which  the  samples 
evolved  subject  to  constant  temperature  conditions  imposed 
using  the  Groningen  method. 
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Spectral  Calculations 


From  the  discussion  in  Chapter  2,  the  expressions  for  the 
isotropic  and  anisotropic  (depolarized)  intensities  are: 


These  quantities  were  calculated  from  the  configurations 
generated  by  the  molecular  dynamics  simulations  described  in 
the  preceding  section.  This  requires  the  transformation  of 
the  body-fixed  frame  components  of  the  susceptibility  tensors 
into  their  laboratory-frame  components  which  was  performed 
using  the  body-fixed  frame  for  water  in  which  the  principal 
moments  of  inertia  are  ordered  < I22  < I33. 

Three  scalar  quantities  may  be  constructed  from  the 

components  of  a,  specifically^®: 


[2-54] 


[2-53] 


[4-5] 


i.e.,  the  molecular  isotropic  point  polarizability,  and 
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[4-6] 


[4-7] 


which  are  related  to  the  optical  anisotropy  6 by 

52  = 0(2  + 3c(2 


2 

2 


[4-8] 


The  molecular  polarizability  of  the  i"*'  molecule  in  the 
laboratory  frame  may  then  be  expressed  as  the  sum  of  its  trace 
and  two  traceless  tensors,  and  according  to: 


in  which  e^n„  (n  = 1,2,3)  is  the  direction  cosine  of  the  n*"^ 
body-fixed  axis  of  the  i*'*’  molecule  with  respect  to  the  a-axis 
of  the  laboratory  frame. 

Observe  that  this  permits  all  tensors  in  the 

spectr\im  (see  equation  2-52)  to  be  represented  as  vectors 
directed  as  axis  2 in  the  body-fixed  frame,  that  is. 


= oc^I  + oc^Q^  + cx^R^ 


[4-9] 


where 


[4-10] 


[4-11] 
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[4-12] 


Ey  PavY  = P ® 


2a 


1 


[4-13] 


The  contribution  of  the  dipole-quadrupole  polarizability  were 
calculated  from  the  components  of  A in  the  laboratory  frame 
constructed  from 


The  values  of  °\i,  3,  and  A are  collected  in  Table  4-3. 

It  should  be  noted  that  these  simulations  used  the  water 
polarizability  components  calculated  by  Huiszoon  using  ab 
initio  methods^®.  The  resulting  polarizability  tensor  is 
significantly  more  anisotropic  than  the  experimental  values 
obtained  by  Murphy®°  and  used  by  de  Santis  et  al  in  their 
calculations  of  Raman  scattering  from  simulated  water.  The 
more  anisotropic  tensor  was  found  by  Bosma  et  al®^  to  give  a 
much  more  realistic  description  of  the  scattering  from 
librational  modes. 

The  time  evolution  of  the  different  polarizability  terms 
may  be  followed  along  the  molecular  dynamics  trajectory  by 
constructing  the  signals^® 


e e „e  A 

iij^a  ii2P  /I3Y  113112^3 


[4-14] 


[4-17] 
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( t)  = Ei  ( t)'  + 0£2i^4  ( t)  ] 


IDID 


dep 


(t)  = ^E„3  Ei  Ej,i  (^j) 


[4-15] 

[4-16] 


for  the  allowed  rotational  and  IDID  contributions  and 


= l«„Ei  Ei.i  E.„v  te»  f-^J) 


3Tm 


^i.o  ( t ) = -^Ei  Ej,i  E„3y  Pi3  ^0,3  ’ ( t ) 


[4-19] 

[4-18] 


for  the  first  isotropic  contributions.  It  should  be  noted 
that  the  second  order  propagator,  T‘^’,  is  commonly  referred  to 
as  the  dipole  propagator  and  may  be  written  in  the  more 
familiar  form"*^: 


T 


(2) 


(3r  . .r  . . 

ij  ij 


[4-20] 


The  individual  spectra  and  their  cross  terms  are  calculated 
according  to. 


^^J(v)  = FT[  ^S(t)]-FT[  ^5(t)]*  + c.c  [4-21] 


where  a and  b represent  any  of  the  spectral  labels  previously 


defined. 


Ill 


Results 

Several  atomic  pair  distribution  functions  were 
calculated  in  order  to  characterize  structural  features  of  the 
simulated  systems.  Though  no  simulations  of  liquid  phase 
hydrate-forming  systems  have  been  found  in  the  literature, 
published  reports  of  solid  hydrate  simulations  are  available 
for  comparison.  Figure  4-1  reports  the  oxygen-oxygen  host 
lattice  pair  distribution  function,  goo(R)f  for  the  methane 
hydrate.  It  may  be  compared  to  Figure  4-2,  from  the 
simulation  by  Tse  et  al'*®  who  obtained  comparable  results. 
There  follows  Figure  4-3,  and  4-4  which  report  the  host 
lattice  pair  distribution  functions  goH(R)  gHH(R) 
respectively. 

Guest-host  pair  distribution  functions  were  also 
calculated  and  Figure  4-5  and  4-6  report  representative 
results  obtained  from  the  solid  methane  hydrate  simulation  and 
the  liquid  phase  methane/  water  simulation  respectively. 
Specifically,  these  figures  characterize  the  structural 
relationship  between  the  methane  guest' s center  and  the  water 
oxygen  atom.  Analagous  quantities  for  the  ethylene  oxide 
hydrate  and  ethylene  oxide/  water  system  are  reported  in 
Figures  4-7  and  4-8. 
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Calculated  spectra  are  presented  and  discussed  in  Chapter 

5. 


Table  4-3.  Molecular  polarizabilities  in  e.s.u. 


Water 

Value 

Reference 

«xx 

1.468 

X 

10-24 

59 

Ctyy 

Ctzz 

1.528 

X 

CM 

1 

O 

1.415 

X 

10-24 

V 

1.855 

X 

O 

1 

CO 

62 

3.130 

X 

10-31 

63 

-^123,1 

1.164 

X 

10-36 

64 

-^123,3 

-0.36 

X 

10-32 

Methane 

tto 

2.593 

X 

10-24 

65 

V 

0 

Ethylene  Oxide 

^XX 

66 

ctyy 

«zz 

V 

1.88 

X 

10-1® 

67 

8 


6 - 

Cn  4 .. 

2 - 

0 -I ' ■ \ 

0 1 2 3 4 5 6 

R (A) 

Figure  4-1.  goo(R)  calculated  from  methane  hydrate 
simulation,  T = 145  K. 


Figure  4-2.  goo(R)  calculated  from  simulations  by  Tse 
et  al . (a)  ice  Ic  at  69  K,  (b)  empty  hydrate  at  110  K 

and  (c)  methane  hydrate  at  145  K. 
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Figure  4-3.  goH(I^)  calculated  from  methane  hydrate 
simulation. 
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Figure  4-4. 
simulation. 
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Figure  4-5.  Guest-host  atom-atom  distribution 
function,  0[(CH2)20]  — 0(H20),  calculated  from  the 
ethylene  oxide  hydrate  simulation 


Figure  4-6.  Guest-host  atom-atom  distribution 
function,  0[(CH2)20]  — 0(H20),  calculated  from  the 
ethylene  oxide/  water  simulation 
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Figure  4-7.  Guest-host  atom-atom  distribution 
function,  C(CH4)  — 0(H20),  calculated  from  the  methane 
hydrate  simulation 


Figure  4-8.  Guest-host  atom-atom  distribution 
function,  C(CH4)  — 0(H20),  calculated  from  the 
ethylene  oxide/  water  simulation 


CHAPTER  5 

ANALYSIS  AND  DISCUSSION 
The  High  Wavenumber  Shift  Region 
Fitting  Procedures 

All  spectra  were  smoothed  once  with  a Savitsky-Golay 
three-point  smoothing  function  before  analysis.  Fitting  of 
spectra  was  performed  using  software  entitled  "RAMAN."®®  The 
The  code  fits  combinations  of  Lorentzian  and  Gaussian  peaks  to 
Raman  spectra  with  a polynomial  background  using  the 
Levenberg-Marquardt  method.  Figure  5-1  indicates  the  results 
of  using  this  procedure  on  a representative  spectral  region. 
Spectral  Features  due  to  THF 

Tetrahydrofuran  contains  several  vibrational  modes  which 
contribute  to  infrared  or  Raman  spectra,  or  to  both,  and  the 
assignments  made  by  Palm,  et  al®®,  are  collected  in  Table  5-1. 
One  would  expect  attractive  intermolecular  interactions 
between  THF  and  water  to  occur  primarily  between  either  of  the 
hydrogen  atoms  on  water  and  the  oxygen  atom  on  THF.  Indeed, 
those  spectral  features  which  derive  from  THF  vibrational 
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Figure  5-1.  Typical  results  fitting  an  experimental  spectrum  (total  Raman 
scattering  from  H20(l)  at  25  °C)  . 
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5-1.  Raman  and  Infrared  Spectral  Data  and 

Assignments®® . 
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modes  primarily  involving  its  aliphatic  portions  are 
essentially  insensitive  to  the  variations  in  temperature  and 
composition  that  are  reported  here.  In  contrast,  vibrational 
motions  of  the  THF  molecule  which  involve  the  stretching  of  C- 
0-C  bonds  give  rise  to  spectral  features  that  are  quite 
sensitive  to  changes  in  the  system. 

Specifically,  the  peak  which  occurs  at  913  cm"^  in  the 
pure  liquid  THF  spectrum  at  25  °C  can  be  attributed  to  ring 
stretching  modes,  denoted  COC  stretch  in  this  discussion. 
With  the  addition  of  water  a shoulder  appears  on  the  low 
wavenumber  side  of  the  COC  stretch  peak.  These  two  features 
are  particularly  important. 

Figure  5-2  illustrates  the  effect  of  system  composition 
on  the  position  of  these  peaks:  the  dilution  of  THF  with  water 
shifts  the  shoulder  to  significantly  lower  wavenumbers.  A 
similar  trend,  but  with  smaller  magnitudes,  is  found  when  the 
weaker  hydrogen  bonding  solvents  methanol  and  chloroform  are 
used. 

Changes  in  system  composition  also  alter  the  ratio  of 
integrated  intensities  for  these  two  THF  peaks  as  is 
illustrated  in  Figure  5-3.  The  effect  of  temperature  on  peak 
position  is  much  weaker  as  can  be  seen  in  Figure  5-4. 
Similarly,  the  effect  of  temperature  on  the  ratio  of 
integrated  intensities  is  given  in  Figure  5-5. 
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Figure  5-2.  Composition  dependence  of  COC/Shoulder 
peak  positions 
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Figure  5-3.  Composition  dependence  of  COC/Shoulder 
Intensity  ratio. 
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Figure  5-4.  Temperature  dependence  of  COC  and  Shoulder 
peak  positions  (0.056  mole  fraction  THF) . 
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Figure  5-5.  Temperature  dependence  of  integrated 
intensity  ratio  (0.056  mole  fraction  THF) . 
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If  vibrational  and  rotational  motion  are  uncorrelated, 
one  can  determine  the  rotational  correlation  time  by  comparing 
the  widths  of  a particular  peak  in  the  isotropic  and 
depolarized  spectrum.  That  is,  one  can  calculate. 
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-1 

rot 


2nc  ( W 


depol 


w.  } 

ISO 


[5-1] 


where  W refers  to  the  full-width  at  half-maximum  (FWHM)  for 
the  specified  spectrum'^®.  Calculations  of  this  type  were 
performed  for  the  THF  COC  stretch  peak  and  the  results  are 
presented  in  Table  5-2  from  which  it  is  readily  apparent  that 
the  assumption  of  vibrational-rotational  statistical 
independence  is  not  justified. 


Table  5-2.  THF  rotational  relaxation  times  calculated  from 
depolarization  ratio  of  COC  peak. 
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Spectral  Features  due  to  Water 

The  high  wavenumber  spectrum  of  liquid  water  contains 
features  due  to  the  three  normal  vibrational  modes  of  the 
monomer,  all  of  which  are  strongly  modified  by  the  influence 
of  intermolecular  interactions. 

The  intense  and  broad  band  from  2800  cm"^  to  3800  cm'^  is 
due  to  intramolecular  OH  vibrations.  These  vibrations  are 
extremely  sensitive  to  the  presence  and  nature  of  local 
hydrogen  bonding  and  thus  contain  information  concerning  the 
intermolecular  structure  of  the  liquid  phase. 

The  complexity  of  liquid  water  has  given  rise  to  two 
distinct  classes  of  models.  Continuum  models  assume  that  all 
water  molecules  are  four-coordinated  as  in  ice,  but  with  a 
continuous  distribution  of  hydrogen  bond  strengths  and 
configurations^^ . 

Mixture  models  postulate  the  existence  of  a small  niomber 
(nine  or  fewer)  of  species  (e.g.,  dimers,  trimers,  pentamers) 
that  differ  in  the  number  and  configuration  of  hydrogen  bonds. 
The  spectra  of  the  liquid  are  usually  separated  into 
components  which  are  then  attributed  to  different  species. 
However,  the  nixmber  of  species  involved  is  an  unknown 
parameter  which  leads  to  controversial  multi-peak  band 
fitting.  Nevertheless,  this  is  the  approach  that  has 
consistently  garnered  experimental  support. 
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The  analysis  of  the  OH  stretching  region  that  is  offered 
here  closely  follows  that  described  by  Walrafen^^  and  belongs 
in  the  class  of  mixture  models. 

In  the  analysis  offered  by  Walrafen,  the  OH  stretching 
region  is  decomposed  into  components  classified  according  to 
whether  they  lie  above  or  below  an  isobestic  point  at  3425  cm' 
Those  lying  above  this  point  are  attributed  to  hydrogen 
bonded  molecules  while  those  below  are  assigned  to  molecules 
not  hydrogen  bonded.  The  integrated  intensities  may  then  be 
used  to  infer  the  relative  number  of  molecules  which  are 
hydrogen  bonded  to  other  water  molecules  compared  to  those 
that  are  not. 

Following  Walrafen,  the  OH  stretch  region  of  the  spectra 
was  first  corrected  using  the  Bose-Einstein  thermal  population 
factor;  it  was  then  decomposed  into  three  Gaussian  components. 
Figure  5-1  indicates  the  results  of  such  a decomposition. 
Components  centered  at  positions  higher  than  3425  cm'^  were 
assigned  to  hydrogen  bonded  monomers  and  their  integrated 
intensities  combined  to  form  the  numerator  of  a ratio  denoted 
f.  The  combined  integrated  intensities  of  components  centered 
below  3245  cm'^  form  the  denominator  of  f.  This  ratio  is 
interpreted  as  the  mole  fraction  of  partially  covalent 
hydrogen  bonded  monomers.  In  Figure  5-7  these  results  are 
plotted  as  a function  of  temperature  where  they  are  compared 
to  those  obtained  by  Walrafen. 
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Walrafen 
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Figure  5-7.  Temperature  dependence  of  hydrogen  bonding 
in  water  compared  to  THF  solution 
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Figure  5-8.  Composition  dependence  of  hydrogen 
bonding  at  25  °C. 
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Aqueous  solutions  of  THF  contain  additional  spectral 
features  that  arise  from  CH  stretching  modes  and  which  appear 
in  this  same  wavevimber  region;  thus,  two  additional  Gaussian 
peaks  must  be  included  in  the  fitting  procedure. 

The  Low  Wavenumber  Shift  Region 

The  Raman  spectra  results  begin  with  those  for  the 
simulations  of  clusters.  Figure  5-10  reports  the  total  Raman 
spectra  as  a function  of  cluster  size,  where  they  are  plotted 
in  the  reduced  form, 

R(co)  = (cOq  - co)  co{  1 - exp  [ -hco/2nJcr]  } J®'^(co)  [2-60] 

Figure  5-11  reports  the  variation  of  Raman  spectra  for 
the  pentamer  as  a function  of  temperature. 

The  spectra  of  the  clusters  demonstrate  how  these  Raman 
signals  converge,  as  a function  of  size,  to  that  of  the  bulk. 
As  one  might  expect,  the  dimer  is  the  main  exception  in  the 
series,  where  the  lowest-  and  highest-f requency  peaks  are 
redshifted  relative  to  the  pentamer  spectrum.  This  is  in 
contrast  to  the  overall  trend  for  redshifting  of  these  peaks 
with  increasing  size.  This  behavior  may  be  the  result  of  the 
two  molecules  in  the  dimer  having  a larger  range  of  rotational 
mobility. 


Reduced  Intensity 
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Figure  5-9.  Reduced  Raman  spectra  calculated  for 
clusters  of  various  sizes. 


Reduced  Intensity 
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Figure  5-10.  Reduced  Raman  spectra  calculated  fron 
water  pentamer  simulations  at  various  temperatures. 
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It  is  doubtful  that  the  peak  which  appears  at  45  cm"^  in 
the  pentamer  has  the  same  origin  as  the  60  cm"^  peak  in  the 
bulk  since  the  former  redshifts,  as  a function  of  size  from  45 
cm"^  in  the  pentamer  spectrum  to  25  cm”^  in  the  simulation  of 
50  molecules.  As  size  is  increased,  the  peaks  which  are  near 
4 5 cm"^  and  700  cm'^  in  the  pentamer  either  diminish  or 
broaden,  and  shift  to  the  red.  The  700  cm'^  peak  is  not 
clearly  resolved  in  the  larger  clusters.  The  200  cm'^  peak  has 
a distinct  bimodal  structure  in  the  smaller  clusters 
suggesting  that  more  than  one  type  of  motion  may  be 
contributing  to  this  feature. 

There  are  several  features  of  note  in  the  spectra  of  the 
pentamer  simulations,  presented  as  a function  of  temperature 
in  Figure  5-11.  The  peak  near  45  cm'^  has  a distinct  bimodal 
character  at  lower  temperatures.  The  peak  broadens  only 
slightly  as  the  temperature  increases  and  it  shifts  to  the  red 
only  a few  wavenumbers. 

The  peak  near  140  cm”^  broadens  considerably  with 
temperature  increase.  A shoulder  on  its  high  frequency  side 
(150-170  cm"^)  increases  in  intensity  with  temperature 
increase,  surpassing  the  140  cm'^  peak  near  room  temperature. 

The  librational  features  can  be  seen  to  undergo 
significant  redshifts  as  the  temperature  increases.  This  can 
be  understood  in  that  the  molecules  have  more  energy  to  move 
at  the  higher  temperatures,  and  consequently  have  larger 
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amplitude  librational  motions  with  lower  characteristic 
oscillation  frequencies. 

Raman  spectra  for  the  bulk  water  simulations  are 
decomposed  into  their  components  and  left  in  their  unreduced 
form  begining  with  Figure  5-12  which  contains  the  allowed 
rotational  spectra,  reported  as  a function  of  temperature. 
These  show  the  characteristic  lineshapes  of  the  orientational 
contributions  of  asymmetric  tops,  specifically,  an  intense  Q- 
branch  plus  a broad  band  extending  past  700  cm'^.  They  also 
show  a well-defined  band  arround  480  cm"^  in  the  supercooled 
liquid.  With  an  increase  in  temperature  this  band  is  red- 
shifted  until  it  merges  into  the  Q-branch. 

From  Figure  5-13  it  can  be  seen  that  the  allowed  portion 
of  the  depolarized  spectrum  is  dwarfed  by  the  IDID 
contribution  which  is  at  least  an  order  of  magnitude  larger. 

Since  tto  is  a scalar,  the  dynamical  information  appears  in  the 

IDID  spectrum  through  the  dipole  propagator  and  the  spectra 
reflect  only  the  center  of  mass  translational  dynamics.  In 
simple  fluids  the  lineshapes  are  approximately  exponential, 
while  water  spectra  develop  more  structured  profiles 
containing  features  which  are  probably  related  to  those 
observed  in  inelastic  neutron  experiments  and  in  the  power 
spectrum  of  the  center  of  mass  velocity. 
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Figure  5-13.  Higher  order  contributions  to  the  isotropic 
spectrum  calculated  for  the  288  K bulk  watersimulation. 


Figure  5-14.  The  total  isotropic  spectrum  calculated 
from  the  288  K bulk  water  simulation. 
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Figure  5-15.  Representative  depolarized  experimental 
water  spectra. 


Figure  5-16.  Representative  depolarized  spectra 
obtained  from  THF/  water  solutions  (0.056  mole 
fraction  THF) . 
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Figure  5-17.  Depolarized  spectra  calculated  from  bulk 
water  simulations. 


Reduced  Intensity 
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Figure  5-18.  Depolarized  spectra  calculated  from 
ethylene  oxide/water  simulations. 


Reduced  Intensity 
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Figure  5-19.  Depolarized  spectra  calculated  from 
ethylene  oxide/water  simulations. 


141 


Higher  order  contributions  to  the  light  scattering 
spectra  are  plotted  for  the  288  K MD  simulation  in  Figure  5- 
14.  Though  these  are  calculated  from  different  tensors,  all 
of  them  show  an  approximately  exponential  decay  to  200  cm‘^ 
with  a characteristic  frequency  of  about  80  cm'^.  This  is 
followed  by  a broad  tail  which  extends  to  about  900  cm‘^. 

The  sum  of  these  three  higher  order  contributions  is 
plotted  in  Figure  5-15  as  the  total  isotropic  spectrum.  Its 
lineshape  appears  to  be  independent  of  temperature.  Like 
their  experimental  counterparts,  the  calculated  isotropic 
spectra  are  essentially  featureless.  The  calculated  spectra 
suffer  an  additional  handicap:  their  dependence  on  the  poorly 
known  dipole-quadrupole  polarizability,  A.  This  quantity  is 
known  for  only  a handful  of  molecules,  and  no  experimental 
values  for  the  dipole-quadrupole  polarizability  could  be  found 
for  either  methane  or  ethylene  oxide;  calculated  values  range 

over  nearly  two  orders  of  magnitude.  Since  the  aTA  terms  are 

of  the  same  order  as  the  other  higher-order  contributions, 
their  neglect  would  lead  to  unreliable  isotropic  spectra; 
hence  the  remainder  of  this  discussion  will  concern  only  the 
depolarized  spectra. 

Both  the  experimental  and  simulated  low  wavenumber 
depolarized  spectra  for  bulk  water  are  consistent  with  the 
results  obtained  by  others.  Specifically,  they  contain  maxima 
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located  at  nearly  the  same  positions,  and  show  the  correct 
temperature  dependence. 

Decomposition  of  the  room  temperature  experimental 
reduced  depolarized  spectrum  into  Gaussian  components  gives 
components  centered  at  approximately  65  cm"^  and  190  cm"^, 
assigned  to  restricted  translational  modes  of  motion,  and  very 
broad  components  centered  near  470,  550,  and  770  cm"^,  which 
are  assigned  to  librational  motion.  The  simulated  room 
temperature  spectrum  yields  nearly  identical  results:  65,  195, 
477,  609,  and  752  cm‘^.  These  compare  favorably  to  the 
experimental  Raman  results  reported  by  Walrafen^^:  60,  170, 
470,  550,  775  cm'^ . In  the  spectrum  of  water  far-infrared 
absorption  coefficients^"*,  broad  bands  were  found  at  60,  and 
188  cm"^  with  an  unresolved  librational  maximum  at 
approximately  685  cm"^.  Similar  results  have  been  obtained  in 
calculations  of  the  vibrational  density  of  states  derived  from 
inelastic  neutron  scatterings^. 

As  Figure  5-17  makes  clear,  the  addition  of  THF  produces 
rather  large  effects,  both  in  the  position  and  intensity  of 
the  maxima,  particularly  in  the  restricted  translation  region. 

Figure  5-18  contains  the  reduced  depolarized  spectra 
obtained  from  bulk  water  simulations.  An  effect  similar  to 
that  produced  by  the  addition  of  THF  in  the  experimental 
spectra  also  appears  in  the  spectra  calculated  from  the 
simulations,  the  ethylene  oxide/  water  system,  given  in  Figure 
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5-19.  This  is  followed  by  the  simulated  methane/  water  system 
spectra  given  in  Figure  5-20. 

Compared  to  their  high  wavenumber  counterparts,  the  low 
wavenumber  spectra  are  characterized  by  broader  and  much  less 
structured  features.  This  makes  their  interpretation  more 
difficult  but  some  insight  can  be  obtained  from  power  spectra 
calculated  from  the  simulated  solid  hydrates. 

The  velocity  time  autocorrelation  function  is  given  by. 


C (t) 

V ' ' 


<v(  t)  • v(  0)  > 
<v(0)  • v(0)  > 


[5-2] 


and  the  power  spectrum  associated  with  translational  motion  is 
obtained  as  the  Fourier  transform  of  the  velocity 
autocorrelation  function^®.  For  a monatomic  harmonic  crystal, 
the  power  spectrum  is  just  the  density  of  vibrational  states. 
For  a molecular  system,  one  may  study  analogous  quantities 
describing  the  translational  motion  of  the  center  of  mass  and 
rotational  motion  of  the  molecule,  i.e.,  the  autocorrelation 
functions  of  the  molecular  velocity  and  angular  velocity, 
respectively. 

A good  approximation  of  the  host  lattice  translational 
density  of  states  is  given  by  the  power  spectrum  of  oxygen 
atom  motions  and  this  is  presented  in  Figure  5-21  for  the 
simulated  ethylene  oxide  hydrate.  Similarly,  one  may  obtain 
insight  into  the  host  librational  motion  from  Figure  5-22 
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which  displays  the  power  spectriom  obtained  from  the  motion  of 
hydrogen  atoms  of  the  host  lattice.  Analogous  spectra 
corresponding  to  motion  of  the  guest  are  presented  in  Figure 
5-17  and  5-18. 

These  calculations  were  also  performed  for  the  simulated 
methane  hydrate.  The  motion  of  host  water  molecules  is 
similar  in  nature  in  the  two  hydrates  though,  as  Table  5-2 
indicates,  there  are  differences  in  the  details.  The  motion 
of  the  guest  molecules  in  the  two  hydrate  types  differs  more 
substantially,  as  Figures  5-23  through  5-26  make  clear. 

The  main  features  of  the  translational  density  of  states 
for  the  hydrates  are  similar  to  those  found  in  experimentally 
for  ice  Ic  at  a comparable  temperature^^.  Even  for  the  solid 
hydrate,  variations  in  temperature  lead  to  a significant 
redistribution  of  the  density  of  states.  Specifically,  an 
increase  in  temperature  leads  to  a very  large  enhancement  in 
the  density  of  states  in  the  100-300  cm"^  region  of  the 
spectra.  Furthermore,  there  is  a tendency  for  the  vibrations 
to  shift  to  lower  temperatures.  This  trend  has  been  observed 
in  the  infrared  spectra  of  ice  Ic  and  ice  Ih^®,  in  the  infrared 
and  Raman  spectra  of  liquid  water^®,  and  in  the  bulk  liquid 
simulations  reported  here. 

An  important  feature  of  the  translational  spectrum  is  the 
well  defined  gap  appearing  at  approximately  116  cm‘^  and  140 
cm"^  for  the  methane  and  ethylene  oxide  hydrates  respectively. 
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This  gap,  though  not  as  well  defined,  appears  in  the 
vibrational  power  spectrum,  the  Raman  spectrum  of  the 
simulated  liquid  systems  and  in  the  experimental  THF/  water 
Raman  spectrum,  though  at  significantly  higher  wavenumbers. 
These  may  be  compared  to  the  results  of  incoherent  inelastic 
neutron  scattering  for  the  SI  SO2  hydrate®®,  whose  low 
frequency  neutron  time-of-f light  spectrum  showed  two  broad 
maxima  at  69  and  244  cm'^,  separated  by  dip  at  110  cm"^.  The 
large  scattering  cross  section  of  the  proton  prevented 
observation  of  SO2  motion. 

The  translational  frequency  spectrum  of  the  encaged 
methane  molecules  contains  well  defined  peaks  at  35,  54,  and 
78  cm"^.  By  following  the  motion  of  specific  molecules,  these 
can  be  partitioned  into  the  contribution  from  the  small  cages 
and  large  cages,  mean  radii  ~4.0  and  ~4.3  A respectively. 
Thus  the  higher  energy  vibration,  78  cm"^,  is  assigned  to  the 
motion  of  methane  molecules  in  the  smaller,  nearly  spherical 
cages.  The  lower  vibrational  frequencies  and  the  splitting  of 
the  translational  energy  of  guests  in  the  large  cages  may  be 
understood  as  a consequence  of  the  greater  size  and  assymmetry 
of  the  cavity.  The  fact  that  these  frequencies  fall  in  the 
same  range  as  several  prominent  peaks  in  the  host 
translational  spectra  suggests  that  even  in  the  solid  hydrate 
the  potential  exists  for  coupling  between  the  translational 
motion  of  the  guest  molecules  and  that  of  the  host  lattice. 
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Although  there  are  no  experimental  infrared  data 
available  for  methane  hydrates,  there  are  results  for  the 
xenon  hydrate  between  50  and  350  cm"^®^.  The  main  features  are 
similar  to  those  for  ice  !<,,  though  there  is  a slight  shift  to 
higher  energy  in  the  hydrate,  which  is  in  general  accord  with 
the  results  reported  here. 

With  respect  to  the  librational  motion  of  the  host 
lattice,  the  calculated  frequency  spectrum  spans  the  range 
between  450  and  1000  cm"^.  Proton  motion  can  be  studied 
experimentally  by  incoherent  neutron  scattering,  and  for  ice 
Ih,  the  neutron  time-of-f light  cross  section  shows  a broad 
asymmetric  profile  with  a maximum  located  at  57  0 cm"^  and  a 
total  width  of  approximately  600  cm'^  ®^. 

In  contrast  to  the  pattern  found  for  the  host  lattice, 
the  librational  modes  of  the  guest  occur  at  lower  frequencies 
than  the  translations.  The  methane  molecules  trapped  in  the 
clathrate  rotate  essentially  freely  about  all  inertial  axes. 
The  librational  density  of  states  for  ethylene  oxide  molecules 
is  more  complicated,  perhaps  due  to  strong  coupling  with  its 
translational  motion.  Two  absorption  peaks  located  at  65.0 
and  35.6  cm"^  have  been  observed  in  the  infrared  spectrum  of 
the  ethylene  oxide  hydrate®®,  suggesting  that  the  four  peaks 
appearing  in  the  librational  density  of  states  should  be 
interpreted  as  a single  pair  which  are  split  into  doublets  due 
to  long  range  dipole-dipole  interactions.  If  this  is  the 
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case,  then  the  infrared  features  should  be  compared  to  the 
mean  values  of  61  and  27  cm"^  calculated  from  the  simulation. 

A neutron  inelastic  scattering  experiment®^  yielded  three 
main  peaks  at  29,  39,  and  57  cm'^.  Good  agreement  with 
dielectric  and  NMR  results  for  the  activation  energy  was 
obtained  if  the  lowest  frequency  was  assigned  to  libration  of 
the  (CH2)20  molecule  about  its  polar  axis  and  the  higher 
frequency  librations  to  motions  about  the  inertia  axes  which 
are  perpendicular  and  parallel  to  the  dipole  axis, 
respectively. 

To  discuss  the  reorientational  motion  of  the  molecules 
one  may  calculate  the  time  autocorrelation  functions  Ci*(t)  of 
the  first  four  Legendre  polynomials  (1=1,  2,  3,  4)  of  the 
angle  subtended  by  the  molecular  axis  corresponding  to  the 
dipole  axis  with  respect  to  its  position  at  t=0. 

C/(t)  = <P^[e*(t)-e^(0)  ]>  [5-3] 

In  Figure  5-26  the  results  from  the  288  K bulk  water 
simulations  are  plotted  in  the  form  [1(1  + 1)]”^  In  Ci*(t). 
This  type  of  presentation  has  been  found  useful  in  evaluating 
the  merits  of  models  of  molecular  reorientation. 

The  four  functions  plotted  in  Figure  5-26  have  the  same 
basic  features.  Specifically,  a glitch  appears  at  times  so 
short  that  they  would  correspond  to  almost  free  inertial 
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motion  in  simple  molecular  fluids.  This  is  followed  by  the 
onset  of  the  diffusive,  slowly  decaying  exponential  regime. 

Though  the  motion  at  long  times  may  be  described  as 
diffusive,  the  behavior  is  not  Debye-like  since  the  functions, 
as  plotted,  do  not  coincide.  The  form  in  which  the  results 
are  shown  also  disguises  the  fact  that  the  amplitude  of  the 
glitch  increases  rapidly  with  increasing  1. 

In  most  dynamical  models  of  liquid  water,  each  molecule 
is  assumed  to  reside  as  part  of  a relatively  long-lived  cage 
of  four  tetrahedrally  coordinated  neighbors.  This  is  the  same 
approach  used  by  Walrafen  to  describe  the  low  frequency  Raman 
scattering  of  liquid  water.  In  this  approach  the  glitch  in 
the  orientational  correlation  functions  is  taken  to  reflect 
the  librational  motion  of  the  central  molecule  in  the  cage  and 
can  therefore  be  linked  to  the  prominent  librational  features 
which  appear  in  both  the  infrared  and  Raman  spectra.  In 
contrast,  the  exponential  regime  is  associated  with  the 
structural  breakup  of  the  cage.  The  rate  of  decay  largely 
determines  the  NMR  relaxation  time  and  its  smooth  character 
accounts  for  the  simplicity  of  the  dielectric  dispersion  in 
the  microwave  region®®. 

In  Figure  5-27  the  1=2  autocorrelation  functions  for  the 
bulk  water  simulations  at  288  K are  plotted  for  the  molecule- 
fixed  X,  y and  z directions.  The  three  functions  have  the 
same  general  shape  and  the  position  of  the  glitch  is  nearly 
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the  same  in  every  case,  in  fact  the  y and  z functions  are 
nearly  identical. 

Orientational  correlation  functions  such  as  these  were 
fit  to  an  exponential  at  large  t (0.8  to  1.5  ps)  from  which  an 
estimate  of  the  reorientational  relaxation  time  is  obtained. 
These  have  been  tabulated  (Table  5-3)  along  with  calculated 
self-diffusion  coefficients  which  were  obtained  by  monitoring 
the  displacement  of  a labelled  atom  i as  a function  of  time, 
averaging  over  all  atoms  and  all  choices  of  time  origin,  and 
exploiting  the  Einstein  relation 

6D  = lim,  _^<|i?.(t)  -i?.(0)|2>  [5-4] 

The  agreement  with  experimental  data  is  only  fair,  the 
calculated  reorientation  times  are  consistently  high  but  their 
temperature  dependence  is  relatively  good. 

The  peaks  which  appear  in  the  calculated  Raman  spectra  at 
frequencies  greater  than  about  300  cm"^  arise  from  the  same 
physical  processes  which  lead  to  the  glitch  in  the 
orientational  correlation  function.  These  features  are  in 
good  agreement  with  experiment  in  several  respects.  First, 
the  peak  positions  are  almost  independent  of  temperature. 
Experimentally,  this  observation  has  led  to  the  suggestion 
that  liquid  water  contains  clusters  of  well-defined  structure. 
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However,  they  may  be  associated  with  short  lived  features  in 
the  orientational  correlation  function. 

Secondly,  the  decrease  in  intensity  of  the  bands  with 
temperature  increase  is  in  semi-quantitative  agreement  with 
experiment.  This  reflects  the  fact  that  the  depth  of  the 
glitch  decreases  as  the  temperature  rises  because  the  cage 
surrounding  a given  molecule  relaxes  more  rapidly  at  high 
temperatures . 

The  suggestion  that  the  low  frequence  Raman  scattering  of 
liquid  may  indicate  well  defined  water  clusters  has  been 
developed  by  Walrafen®^.  Figure  5-28,  modified  from  his 
discussion,  displays  a five-molecule,  fully  hydrogen-bonded 
structure  taken  from  the  ice  Ij,  lattice.  The  overall  symmetry 
of  the  structure  is  C2.  However,  in  terms  of  the  central 
oxygen  atoms,  the  tetrahedral  nearest-neighbor  oxygen  atoms, 
and  the  four  intervening  protons,  the  symmetry  is  C2v»  Then 
the  symmetry  may  be  described  as  intermolecular  C2V  symmetry, 
provided  the  conformations  of  the  outer  protons  are  excluded. 

In  this  model  nine  intermolecular  vibrations  must  be 
considered,  the  intermolecular  vibrations  being  excluded.  Of 
these,  six  refer  to  restricted  translations  and  three  to 
restricted  rotations,  i.e.,  librations. 

Two  of  the  six  restricted  translations  are  equivalent  to 
the  hydrogen  bond  bending  motions.  They  are  the  totally 
symmetric  deformations  belonging  to  the  vibrational  species 
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and  are  designated  VjA^  and  v^Aj . The  remaining  four  restricted 
translations,  hydrogen  bond  stretching  motions,  involve  the 
Ai,  Bj,  and  B2  vibrational  species  and  are  designated  v^Ai,  V2A1, 

V7B1,  and  V9B2.  The  hydrogen  bond  bending  vibrations  V3A1  and 

V4A1  are  unresolved  and  occur  near  60  cm”^;  the  hydrogen  bond 

stretching  vibrations  VjAi,  V2A1,  V7B1,  and  V9B2  are  also 

unresolved  and  give  rise  to  the  broad  band  centered  near  170 
cm"^  for  liquid  H2O. 

The  restricted  rotational  components  include  one 
designated  V5A2  which  involves  the  restricted  rotation  that  is 
symmetric  only  with  respect  to  the  C2  axis.  It  is  forbidden 
in  all  dipolar  spectra,  namely  infrared  and  Pi  hyper-Raman 

spectra.  The  band  designated  VgB2  is  assigned  to  the  libration 

that  is  symmetric  with  respect  to  the  plane  formed  by  the 
atoms  of  the  central  H2O  molecule  of  figure  5-28.  That  plane 
is  arbitrarily  designated  B2 . The  remaining  libration, 

designated  VgBj  is  symmetric  only  with  respect  to  the  Bi  plane 

(perpendicular  to  the  B2  plane  and  passing  through  the  C2 
axis) . Because  none  of  the  librations  is  totally  symmetric  to 
all  symmetric  elements  (C2,  Bi,  B2)  , each  is  completely 
depolarized  in  the  Raman  spectrum  for  C2V  symmetry. 
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Figure  5-20.  Power  spectrum  of  water  oxygen  atoms 
obtained  from  simulated  ethylene  oxide  hydrate. 
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Figure  5-21.  Power  spectrum  of  water  hydrogen  atoms 
obtained  from  simulated  ethylene  oxide  hydrate. 
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Figure  5-22.  Power  spectrum  of  ethylene  oxide 
translational  motion  obtained  from  simulated 
hydrate . 
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Figure  5-23.  Power  spectrum  of  ethylene  oxide 
libration  obtained  from  simulated  hydrate. 
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Figure  5-24.  Power  spectrum  of  methane  translation 
obtained  from  simulated  hydrate. 


Figure  5-25.  Power  spectrum  of  methane  libration 
obtained  from  simulated  hydrate. 
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For  liquid  H2O  the  V5A2  librations  are  centered  near  450 
the  VgBi  librations  are  centered  near  750  cm"^,  and  the 

V8B2  librations  are  centered  near  550  cm"^. 

The  structure  is  useful  in  making  vibrational 
assignments,  but  it  necessarily  refers  to  an  average 
configuration.  Distortions  from  C2V  symmetry  involving  angles 
and  distances  are  present  in  those  molecules  of  the  real 
liquid  that  are  involved  in  the  fully  hydrogen-bonded 
structures,  and  these  various  distortions  give  rise  to  the 
breadth  observed  for  the  individual  intermolecular  components . 
Distortions  are  also  present  in  the  partially  hydrogen-bonded 
structures . 

A second,  and  even  more  significant,  weakness  in  this 
approach  is  that  it  cannot  be  readily  extended  to  solutions 
containing  components  in  addition  to  water.  This  is 
particularly  true  for  the  systems  considered  here  since  even 
in  the  solid  phase  they  contain  unit  cells  that  are  rather 
large . 

Quantitative  analysis  of  the  experimental  low  frequency 
Raman  spectra  is  made  problematic  by  the  high  noise  levels  and 
the  broad  background  due  to  Rayleigh  scattering  and  stray 
light.  The  spectra  calculated  from  the  simulations  contain 
broad  unresolved  features  which  may  not  necessarily  have  well- 
defined  Gaussian  or  Lorentzian  profiles. 
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However,  it  is  possible  to  take  advantage  of  the  clear 
break  that  occurs  between  the  water  translational  and 
librational  regions  of  the  spectra.  Thus,  the  following 
procedure  was  adopted:  calculated  and  experimental  spectra 
were  decomposed  into  the  minimum  number  of  Gaussian  components 

required  to  achieve  a fit  having  a value  comparable  to  that 

obtained  for  the  experimental  water  spectra;  the  "peaks"  were 
then  classified  as  belonging  to  the  water  translational  or 
librational  region  according  to  whether  they  were  centered 
above  or  below  the  break  in  intensity  found  in  each  spectrum. 
The  integrated  intensities  are  reported  as  a percentage  of  the 
total,  excluding  background  contributions,  in  Table  5-3. 

This  classification  is  is  not  meant  to  indicate  the 
physical  origin  of  specific  spectral  components,  since  guest 
librations  generally  appear  in  the  same  spectral  region  as 
water  translations,  for  example.  More  importantly,  the  DID 
model  presupposses  that  the  relevant  polarizability 
fluctuations  are  due  to  motions  involving  several  molecules 
so  it  is  not  physically  meaningful  to  assign  individual 
features  to  specific  types  of  motion  performed  by  specific 


molecules . 
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Recomendations  for  Further  Study 

Agreement  between  the  experimental  and  computational 
results  is  semi-quantitative.  This  should  not  be  taken  for 
granted  since  the  temperatures  are  low  and  the  systems  include 
light  atoms;  so  that  quantum  effects  should  be  non-negligible. 
There  are  a number  of  ways  in  which  the  experimental  and 
computational  aspects  of  these  studies  can  be  improved  and 
extended. 

The  low  water  soliibility  of  most  clathrate  hydrate  guests 
presents  a challenge  to  the  experimentalists.  Clathrate 
formation  typically  begins  at  the  gas-liquid  interface,  but 
that  interface  is  much  more  efficient  at  scattering  light  than 
the  microscopic  hydrate  nuclei  and  thus  is  an  important  source 
of  noise  in  the  experimental  spectra.  Overcoming  this 
difficulty  would  require  a redesign  of  the  scattering  cell. 
Similarly,  Raman  studies  of  solids  are  best  performed  with  a 
cell  that  rotates  rapidly  such  that  mesoscopic  anisotropies 
are  averaged  and  this  would  also  require  a redesign  of  the 
scattering  cell. 

The  experimental  studies  reported  here  were  obtained 
using  a traditional  monochrometer  that  scans  through  the 
wavelength  range  of  interest.  This  is  a slow  process  and  the 
resulting  spectra  must  be  considered  those  of  systems  at 
equilibrium.  A study  directly  addressing  the  kinetics  of 
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hydrate  formation  would  require  use  of  the  much  faster  Fourier 
transform  method  of  Raman  spectroscopy.  Fourier  transform 
Raman  instruments  also  typically  yield  spectra  of  higher 
resolution  than  available  from  the  apparatus  used  here. 

The  SPC  force  field  used  to  model  water  is  as  good  or 
better  than  any  of  the  other  available  effective-pair 
potentials  for  rigid  water  molecules.  This  is  also  the  case 
for  the  intermolecular  potentials  used  to  model  the  guests. 
Only  marginal  improvement  can  be  expected  from  models  that 
permit  intra-molecular  bond  flexibility  since  the  temperatures 
are  low,  and  the  time-scales  for  intra-molecular  vibrations  in 
the  guests  modeled  here  are  very  fast.  Furthermore,  the  DID 
model  underlying  the  spectral  calculations  is  appropriate  only 
in  the  classical  limit. 

These  models  can  be  significantly  improved  upon  only  by 
incorporating  many-body  or  quantum  effects.  There  is 
currently  a great  deal  of  interest  in  incorporating  many-body 
effects  in  molecular  dynamics  or  Monte  Carlo  simulations  of 
water  and  aqueous  solutions.  Inclusion  of  explicit  three-body 
terms  in  the  potential  function  was  proposed  by  Wojcik  and 
dementi®^,  though  it  is  rather  cumbersome  to  evaluate  and 
parametrize  the  potential.  Early  attempts  to  include  non- 
additive effects  via  a polarizability,  required  the  use  of  a 
time  consuming  interative  procedure  for  determining  the 
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induced  dipole  moments,  which  depend  on  each  other  in  a non- 
linear fashion®®. 

Wallqvist  et  al®®  modified  the  SPC  rigid  water  model  by 
adding  a point  polarizability  and  developed  a new  algorithm 
for  calculating  the  induced  dipole  moments  based  on  a 
predictive,  instead  of  an  iterative  scheme.  This  study 
confirmed  the  importance  of  many-body  effects,  particularly 
with  respect  to  the  description  of  rotational  dynamics,  but 
did  not  significantly  improve  the  description  of  the  structure 
or  internal  energy  of  bulk  water  systems. 

Van  Belle  et  al®°  also  started  with  the  SPC  model  and 
applied  the  extended  Lagrangian  method  to  incorporate  induced 
polarization  effects  in  molecular  dynamics  simulations  of  pure 
water.  The  induced  dipole  on  each  molecule  is  treated  as  an 
additional  degree  of  freedom  and  the  Lagrangian  of  the  system 
is  extended  with  a fictitious  kinetic  term  associated  with 
this  extra  variable.  Their  study  confirmed  a significant 
improvement  in  the  description  of  dynamical  properties  is 
obtained  by  using  a polarizable  model,  and  they  developed  an 
algorithm  that  requires  less  than  twice  the  time  as  what  is 
required  by  an  equivalent  non-polarizable  system. 

Work  on  improving  the  treatment  of  many-body  effects  is 
ongoing,  and,  to  date,  no  polarizable  water  model  has  been 
used  to  calculate  the  Raman  spectrum  of  water  or  any  aqueous 


solution. 
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Figure  5-26.  The  functions  [1(J  + 1)]“^  In  Ci*(t)  for 
the  bulk  water  simulation  at  288  K and  the  values  of  1 
indicated. 


Figure  5-27.  The  functions  C2*(t), 
bulk  water  simulation  at  288  K. 


C2^(t),  C2Mt)  for 
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Figure  5-28.  Five-molecule,  fully  hydrogen-bonded 
structure  having  intermolecular  C2V  symmetry.  Discs 
refer  to  hydrogen  bonds.  Modified  from  reference  73. 
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Table  5-3.  Transport  Properties 

Experimental  Water  Self-Diffusion  Coefficients  and  Rotational 
Relaxation  Times 


T /K 

D /10‘^  s”^  X /ps 

Reference 

259 

8.0 

58 

298 

2.3 

59 

300 

5.0 

60 

350 

0.9 

60 

363 

0.9 

61 

Water 

Times 

Self-Diffusion  Coefficients  and  Rotational 
Calculated  from  Bulk  Water  Simulations 

Relaxation 

T /K 

D /I O'®  s'^ 

X /ps 

255 

3.8 

12 . 8 

288 

4.1 

6.2 

298 

4.3 

5.0 

350 

4.6 

2.0 

Water 

times 

self-diffusion  coefficients  and  rotational  relaxation 
calculated  from  ethylene  oxide/  water  simulations 

T /K 

D /lO'®  s'^ 

X /ps 

290 

3.8 

10.1 

299 

4.1 

6.5 

312 

4.3 

5.4 

330 

4.6 

3.0 

163 


Water  self-diffusion  coefficients  and  rotational  relaxation 
times  calculated  from  methane/  water  simulations 


T /K 

D /lO'®  m®  s"^ 

T /pS 

288 

3.8 

7.3 

300 

4.1 

6.8 

311 

4.3 

4.1 

328 

4.6 

3.5 

Ethylene  oxide 
relaxation  times 

self-diffusion  coefficients  and 
calculated  from  simulations 

rotational 

T /K 

D /I O'®  s'^ 

X /ps 

288 

3.8 

6.0 

300 

4.1 

8.4 

311 

4.3 

7.1 

328 

4.6 

4.3 
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Table  5-4.  Integrated  intensity  below  translational/ 
librational  break  in  depolarized  spectra.  Peported  as 
percentage  of  total  low  frequency  depolarized  integrated 
intensity  excluding  background. 
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